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Abstract 

Energy  resources  and  their  utilization  intimately  relate  to  sustainable  development.  In  attaining 
sustainable  development,  increasing  the  energy  efficiencies  of  processes  utilizing  sustainable  energy 
resources  plays  an  important  role.  The  utilization  of  renewable  energy  offers  a  wide  range  of 
exceptional  benefits.  There  is  also  a  link  between  exergy  and  sustainable  development.  A  sustainable 
energy  system  may  be  regarded  as  a  cost-efficient,  reliable,  and  environmentally  friendly  energy 
system  that  effectively  utilizes  local  resources  and  networks.  Exergy  analysis  has  been  widely  used  in 
the  design,  simulation  and  performance  evaluation  of  energy  systems. 

The  present  study  comprehensively  reviews  exergetic  analysis  and  performance  evaluation  of  a 
wide  range  of  renewable  energy  resources  (RERs)  for  the  first  time  to  the  best  of  the  author’s 
knowledge.  In  this  regard,  general  relations  (i.e.,  energy,  exergy,  entropy  and  exergy  balance 
equations  along  with  exergy  efficiency,  exergetic  improvement  potential  rate  and  some  thermo¬ 
dynamic  parameters,  such  as  fuel  depletion  ratio,  relative  irreversibility,  productivity  lack  and 
exergetic  factor)  used  in  the  analysis  are  presented  first.  Next,  exergetically  analyzed  and  evaluated 
RERs  include  (a)  solar  energy  systems;  (al)  solar  collector  applications  such  as  solar  water  heating 
systems,  solar  space  heating  and  cooling,  solar  refrigeration,  solar  cookers,  industrial  process  heat, 
solar  desalination  systems  and  solar  thermal  power  plants),  (a2)  photovoltaics  (PVs)  and  (a3)  hybrid 
(PV/thermal)  solar  collectors,  (b)  wind  energy  systems,  (c)  geothermal  energy  systems,  (cl)  direct 
utilization  (district  heating,  geothermal  or  ground-source  heat  pumps,  greenhouses  and  drying)  and 
(c2)  indirect  utilization  (geothermal  power  plants),  (d)  biomass,  (e)  other  renewable  energy  systems, 
and  (f)  country  based  RERs.  Studies  conducted  on  these  RERs  are  then  compared  with  the 
previously  ones  in  tabulated  forms,  while  the  Grassmann  (or  exergy  flow)  diagrams,  which  are  a  very 
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useful  representation  of  exergy  flows  and  losses,  for  some  RERs  are  given.  Finally,  the  conclusions 
are  presented.  It  is  expected  that  this  comprehensive  study  will  be  very  beneficial  to  everyone  involved 
or  interested  in  the  exergetic  design,  simulation,  analysis  and  performance  assessment  of  RERs. 

©  2006  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Achieving  solution  to  environmental  problems  that  we  face  today  requires  long-term 
potential  actions  for  sustainable  development.  In  this  regard,  renewable  energy  resources 
(RERs)  appear  to  be  the  one  of  the  most  efficient  and  effective  solutions  [1]. 

RERs  (i.e.,  solar,  hydroelectric,  biomass,  wind,  ocean  and  geothermal  energy)  are 
inexhaustible  and  offer  many  environmental  benefits  compared  to  conventional  energy 
sources.  Each  type  of  renewable  energy  (RE)  also  has  its  own  special  advantages  that  make 
it  uniquely  suited  to  certain  applications.  Almost  none  of  them  release  gaseous  or  liquid 
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Nomenclature 

a  amplitude,  share  of  renewables  in  the  total  of  each  component 

A  area  (m  ) 

C  specific  heat  (kJ/kgK) 

COP  coefficient  of  performance  (dimensionless) 

E  energy  rate  (kW) 

ex  specific  exergy  (kJ/kg) 

Ex  exergy  rate  (kW) 

/  exergetic  factor,  the  sunlight  dilution  factor  (dimensionless) 

F  exergy  rate  of  fuel  (kW) 

g  gravitational  constant  (m/s“) 

h  specific  enthalpy  (kJ/kg) 

HV  heating  value  (kJ/kg) 

HHV  higher  heating  (gross  calorific)  value  (kJ/kg) 

I  global  irradiance  (W/m  ) 

I  rate  of  irreversibility,  rate  of  exergy  consumption  (kW) 

IP  rate  of  improvement  potential  (kW) 

L  enthalpy  of  phase  change  (kJ/kg) 

LHV  lower  heating  (net  calorific)  value  (kJ/kg) 
m  mass  flow  rate  (kg/s) 

n  number,  index  number  (dimensionless) 

P  pressure  (kPa) 

P  exergy  rate  of  the  product  (kW) 

Q  heat  transfer  rate  (kW) 

r  renewable  use  by  the  residential-commercial  sector  in  energy  terms  (kJ) 

R  ideal  gas  constant  (kJ/kgK) 

s  specific  entropy  (kJ/kgK) 

S  entropy  rate  (kW/K) 

SExI  specific  exergy  index  (dimensionless) 

t  period  between  local  maxima  and  minima  of  the  tidal  record,  time  (s) 

T  temperature  (°C  or  K) 

U  heat  transfer  coefficient  (kW/m  K) 

V  speed  (m/s) 

W  work  (kJ) 

W  rate  of  work  (or  power)  (kW) 

y  mol  fraction  (dimensionless) 

z  vertical  distance  from  the  water  level  of  the  reservoir  to  the  reference  height 

or  average  sea  level  (m) 

Z  mass  fraction  (dimensionless) 

Greek  letters 

>y  a 

a  Stefan-Boltzmann  constant  (W/m  K  ) 

k  dilution  factor  (dimensionless) 

rj  energy  (first  law)  efficiency  (dimensionless) 
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flow  (specific)  exergy  (kJ/kg),  maximum  efficiency  ratio  (or  exergy-to-energy 

ratio  for  radiation  (dimensionless) 

interval 

density  (kg/m  ) 

proportionality  constant  (or  quality  factor  or  exergy  coefficient) 

fuel  depletion  rate  (dimensionless) 

exergy  (second  law)  efficiency  (dimensionless) 

productivity  lack  (dimensionless) 

relative  irreversibility  (dimensionless) 

specific  humidity  ratio  (kgwater/kgair) 


air,  actual 

absorber 

adsorbent 

ambient 

atmospheric 

available 

average 

beam 

Carnot 

cooking 

collector-evaporator 

chemical 

collector 

compressor 

condenser 

conversion 

cooker 

cooling 

natural  direct  discharge,  diffuse 
drying  air 

destroyed,  destruction 
diffusive 

electrical,  evaporator 

effective 

engine 

evaporator 

exergetic 

exergetic  residential  commercial 
fuel 

fan-coil 

vaporization 

generator 

geothermal  district  heating  system 
generation,  generated 
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gh 

ground  heat  exchanger 

h 

heating,  heat 

HE 

heat  exchanger 

HHV 

higher  heating  value 

HP 

heat  pump 

HV 

heating  value 

i 

successive  number  of  elements 

ic 

incompressible 

in 

input 

io 

inverter  output 

k 

location 

KN 

kinetic 

LHV 

lower  heating  value 

max 

maximum 

mech 

mechanical 

mix 

mixture 

0 

overall 

oe 

overall  electricity 

of 

overall  fuel 

or 

overall  residential 

out 

output 

P 

constant  pressure,  pump 

par 

parabolic 

per 

perfect 

PH 

physical 

pot 

potential 

PP 

power  plant 

pre 

precipitation 

PT 

potential 

PV 

photovoltaic 

Q 

heat 

r 

reinjected  thermal  water,  refrigerant 

R 

rational 

Ran 

Rankine 

rc 

residential-commercial 

rec 

receiver 

res 

reservoir 

ro 

residential  overall 

s 

solar 

scol 

solar  collector 

sh 

space  heating 

srad 

solar  radiation 

sys 

system 

T 

total,  thermal 

TEG 

thermal  gradient 

trans 

transformation 

u 

useful 
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v  water  vapor 

w  water 

wa  water/antifreeze  solution 

wh  wellhead,  water  heating 

wor  without  reinjected  thermal  water 
wr  with  reinjected  thermal  water 

0  dead  (reference)  state 

1  initial  state 

2  final  state 


pollutants  during  operation.  In  their  technological  development,  the  renewable  ranges 
from  technologies  that  are  well  established  and  mature  to  those  that  need  further  research 
and  development  [1,2]. 

Even  though  conventional  sources,  such  as  oil,  natural  gas  and  coal  meet  most 
of  the  energy  demand  at  the  moment,  the  role  of  RERs  and  their  current  advances  have 
to  take  more  relevance  in  order  to  contribute  to  energy  supply  and  support  the 
energy  conservation  (or  efficiency)  strategy  by  establishing  energy  management  systems 
[3].  The  use  of  RE  offers  a  range  of  exceptional  benefits,  including:  a  decrease  in 
external  energy  dependence;  a  boost  to  local  and  regional  component  manufacturing 
industries;  promotion  of  regional  engineering  and  consultancy  services  specializing  in  the 
utilization  of  RE;  increased  R&D,  decrease  in  impact  of  electricity  production  and 
transformation;  increase  in  the  level  of  services  for  the  rural  population;  creation  of 
employment,  etc.  [4]. 

Dincer  [5]  reported  the  linkages  between  energy  and  exergy,  exergy  and  the  environment, 
energy  and  sustainable  development,  and  energy  policy  making  and  exergy  in  detail.  He 
provided  the  following  key  points  to  highlight  the  importance  of  the  exergy  and  its 
essential  utilization  in  numerous  ways:  (a)  it  is  a  primary  tool  in  best  addressing  the  impact 
of  energy  resource  utilization  on  the  environment,  (b)  It  is  an  effective  method  using  the 
conservation  of  mass  and  conservation  of  energy  principles  together  with  the  second  law  of 
thermodynamics  for  the  design  and  analysis  of  energy  systems,  (c)  It  is  a  suitable  technique 
for  furthering  the  goal  of  more  efficient  energy-resource  use,  for  it  enables  the  locations, 
types,  and  true  magnitudes  of  wastes  and  losses  to  be  determined,  (d)  It  is  an  efficient 
technique  revealing  whether  or  not  and  by  how  much  it  is  possible  to  design  more  efficient 
energy  systems  by  reducing  the  inefficiencies  in  existing  systems,  (e)  It  is  a  key  component 
in  obtaining  a  sustainable  development. 

Sustainable  development  does  not  make  the  world  ‘ready’  for  the  future  generations,  but 
it  establishes  a  basis  on  which  the  future  world  can  be  built.  A  sustainable  energy  system 
may  be  regarded  as  a  cost-efficient,  reliable,  and  environmentally  friendly  energy  system 
that  effectively  utilizes  local  resources  and  networks.  It  is  not  ‘slow  and  inert’  like  a 
conventional  energy  system,  but  it  is  flexible  in  terms  of  new  techno-economic  and  political 
solutions.  The  introduction  of  new  solutions  is  also  actively  promoted  [6]. 

An  exergy  analysis  (or  second  law  analysis)  has  proven  to  be  a  powerful  tool  in  the 
simulation  thermodynamic  analyses  of  energy  systems.  In  other  words,  it  has  been  widely 
used  in  the  design,  simulation  and  performance  evaluation  of  energy  systems.  Exergy 
analysis  method  is  employed  to  detect  and  to  evaluate  quantitatively  the  causes  of  the 
thermodynamic  imperfection  of  the  process  under  consideration.  It  can,  therefore,  indicate 
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the  possibilities  of  thermodynamic  improvement  of  the  process  under  consideration,  but 
only  an  economic  analysis  can  decide  the  expediency  of  a  possible  improvement  [7,8]. 

The  concepts  of  exergy,  available  energy,  and  availability  are  essentially  similar.  The 
concepts  of  exergy  destruction,  exergy  consumption,  irreversibility,  and  lost  work  are  also 
essentially  similar.  Exergy  is  also  a  measure  of  the  maximum  useful  work  that  can  be  done 
by  a  system  interacting  with  an  environment  which  is  at  a  constant  pressure  P0  and  a 
temperature  T0.  The  simplest  case  to  consider  is  that  of  a  reservoir  with  heat  source  of 
infinite  capacity  and  invariable  temperature  T0.  It  has  been  considered  that  maximum 
efficiency  of  heat  withdrawal  from  a  reservoir  that  can  be  converted  into  work  is  the 
Carnot  efficiency  [9,10]. 

Although  numerous  studies  have  been  conducted  on  the  energetic  analysis  and 
performance  evaluation  of  RERs  by  using  energy  analysis  method  in  the  literature,  very 
limited  review  papers  have  appeared  on  exergy  analysis  and  performance  assessment  of 
RERs.  In  this  regard,  Koroneos  et  al.  [11]  dealt  with  the  three  kinds  of  RERs  in  terms  of 
exergetic  aspects,  namely  (i)  exergetic  analysis  of  a  solar  thermal  power  system  by 
presenting  a  study  conducted  by  Singh  et  al.  [12],  (ii)  exergy  analysis  of  geothermal  power 
systems  by  presenting  and  evaluating  the  Larderello-Farinello-Valle  Secolo  Geothermal 
Area  (Tuscany,  Italy)  studied  by  Bettagli  and  Bidini  [13],  and  (iii)  exergy  analysis  of  wind 
energy  systems  by  giving  the  relations  used  in  the  analysis  and  evaluating  a  wind  energy 
system.  They  also  made  a  comparison  between  renewable  and  non-RE  sources,  and 
concluded  that  some  of  the  systems  appear  to  had  high  efficiencies,  and  in  some  cases  they 
are  greater  than  the  efficiency  of  systems  using  non-RE  sources.  In  other  cases,  like  the 
conversion  of  solar  energy  to  electricity,  the  efficiencies  were  lower,  in  order  to  meet  the 
electricity  needs  of  cities.  Hermann  [14]  identified  the  primary  exergy  reservoirs  that  supply 
exergy  to  the  biosphere  and  quantified  the  intensive  and  extensive  exergy  of  their  derivative 
secondary  reservoirs,  or  resources.  Exergy  relations  of  cosmic  radiation  exchange,  wind, 
ocean  surface  waves,  precipitation,  ocean  thermal  gradient,  tides  and  geothermal  were  also 
briefly  presented. 

The  present  study  differs  from  the  previous  ones  due  to  the  facts  that:  (i)  This  covers  a 
comprehensive  exergetic  analysis  and  performance  evaluation  of  RERs.  (ii)  This  includes  a 
wide  range  of  RERs  such  as  solar,  wind,  geothermal  and  biomass  along  with  their 
subsections  (i.e.,  for  geothermal,  its  direct  and  indirect  applications)  as  well  as  hybrid 
systems,  (iii)  This  presents  a  comparison  between  previously  conducted  studies  in  tabulated 
forms,  and  (iv)  This  summarizes  exergetic  utilization  efficiency  of  RERs  as  done  by 
Hepbasli  and  Utlu  [15]  for  the  first  time  for  Turkey.  In  this  regard,  the  structure  of  the 
paper  is  organized  as  follows:  The  first  section  includes  the  introduction;  Section  2  deals 
with  energy  and  exergy  modeling,  while  general  relations  are  described  in  Section  3; 
exergetic  analysis  and  evaluating  the  RERs  are  treated  in  Section  4  in  more  detail  by 
applying  the  general  relations  to  various  RERs;  and  the  last  section  concludes. 

2.  Energy  and  exergy  modeling 

Dincer  et  al.  [16]  reported  that,  to  provide  an  efficient  and  effective  use  of  fuels,  it  is 
essential  to  consider  the  quality  and  quantity  of  the  energy  used  to  achieve  a  given 
objective.  In  this  regard,  the  first  law  of  thermodynamics  deals  with  the  quantity  of  energy 
and  asserts  that  energy  cannot  be  created  or  destroyed,  whereas  the  second  law  of 
thermodynamics  deals  with  the  quality  of  energy,  i.e.,  it  is  concerned  with  the  quality  of 
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energy  to  cause  change,  degradation  of  energy  during  a  process,  entropy  generation  and 
the  lost  opportunities  to  do  work.  More  specifically,  the  first  law  of  thermodynamics  is 
concerned  only  with  the  magnitude  of  energy  with  no  regard  to  its  quality;  on  the  other 
hand,  the  second  law  of  thermodynamics  asserts  that  energy  has  quality  as  well  as 
quantity.  By  quality,  it  means  the  ability  or  work  potential  of  a  certain  energy  source 
having  certain  amount  of  energy  to  cause  change,  i.e.,  the  amount  of  energy  which  can  be 
extracted  as  useful  work  which  is  termed  as  exergy.  First  and  second  law  efficiencies  are 
often  called  energy  and  exergy  efficiencies,  respectively.  It  is  expected  that  exergy 
efficiencies  are  usually  lower  than  the  energy  efficiencies,  because  the  irreversibilities  of  the 
process  destroy  some  of  the  input  exergy. 

Exergy  is  the  expression  for  loss  of  available  energy  due  to  the  creation  of  entropy  in 
irreversible  systems  or  processes.  The  exergy  loss  in  a  system  or  component  is  determined 
by  multiplying  the  absolute  temperature  of  the  surroundings  by  the  entropy  increase. 
Entropy  is  the  ratio  of  the  heat  absorbed  by  a  substance  to  the  absolute  temperature  at 
which  it  was  added.  While  energy  is  conserved,  exergy  is  accumulated  [17]. 

Exergy  analysis  provides  a  method  to  evaluate  the  maximum  work  extractable  from  a 
substance  relative  to  a  reference  state  (i.e.,  dead  state).  This  reference  state  is  arbitrary,  but 
for  terrestrial  energy  conversion  the  concept  of  exergy  is  most  effective  if  it  is  chosen  to 
reflect  the  environment  on  the  surface  of  the  Earth.  The  various  forms  of  exergy  are  due  to 
random  thermal  motion,  kinetic  energy,  potential  energy  associated  with  a  restoring  force, 
or  the  concentration  of  species  relative  to  a  reference  state.  In  order  to  establish  how  much 
work  potential  a  resource  contains,  it  is  necessary  to  compare  it  against  a  state  defined  to 
have  zero  work  potential.  An  equilibrium  environment  which  cannot  undergo  an  energy 
conversion  process  to  produce  work  is  the  technically  correct  candidate  for  a  reference 
state  [14]. 

It  should  be  noticed  that  exergy  is  always  evaluated  with  respect  to  a  reference 
environment  (i.e.  dead  state).  When  a  system  is  in  equilibrium  with  the  environment,  the 
state  of  the  system  is  called  the  dead  state  due  to  the  fact  that  the  exergy  is  zero.  At  the 
dead  state,  the  conditions  of  mechanical,  thermal,  and  chemical  equilibrium  between  the 
system  and  the  environment  are  satisfied:  the  pressure,  temperature,  and  chemical 
potentials  of  the  system  equal  those  of  the  environment,  respectively.  In  addition,  the 
system  has  no  motion  or  elevation  relative  to  coordinates  in  the  environment.  Under  these 
conditions,  there  is  neither  possibility  of  a  spontaneous  change  within  the  system  or  the 
environment  nor  an  interaction  between  them.  The  value  of  exergy  is  zero.  Another  type  of 
equilibrium  between  the  system  and  environment  can  be  identified.  This  is  a  restricted  form 
of  equilibrium,  where  only  the  conditions  of  mechanical  and  thermal  equilibrium  (thermo¬ 
mechanical  equilibrium)  must  be  satisfied.  Such  state  is  called  the  restricted  dead  state.  At 
the  restricted  dead  state,  the  fixed  quantity  of  matter  under  consideration  is  imagined  to  be 
sealed  in  an  envelope  impervious  to  mass  flow,  at  zero  velocity  and  elevation  relative  to 
coordinates  in  the  environment,  and  at  the  temperature  T0  and  pressure  P0  taken  often  as 
25  °C  and  1  atm  [18]. 

3.  General  relations 

For  a  general  steady  state,  steady-flow  process,  the  four  balance  equations  (mass, 
energy,  entropy  and  exergy)  are  applied  to  find  the  work  and  heat  interactions,  the  rate  of 
exergy  decrease,  the  rate  of  irreversibility,  the  energy  and  exergy  efficiencies  [16,19-21]. 
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3.1.  Mass,  energy,  entropy  and  exergy  balances 


The  mass  balance  equation  can  be  expressed  in  the  rate  form  as 


(1) 


where  m  is  the  mass  flow  rate,  and  the  subscript  in  stands  for  inlet  and  out  for  outlet. 

The  general  energy  balance  can  be  expressed  below  as  the  total  energy  inputs  equal  to 
total  energy  outputs. 


y  ^  P'm  — 


(2) 


In  the  absence  of  electricity,  magnetism,  surface  tension  and  nuclear  reaction,  the 
total  exergy  of  a  system  Ex  can  be  divided  into  four  components,  namely  (i)  physical 

•  PH  •  KN  •  PT 

exergy  Ex  ,  (ii)  kinetic  exergy  Ex  ,  (iii)  potential  exergy  Ex  ,  and  (iv)  chemical  exergy 
ExCH  [16]. 


•  PH  •  KN  •  PT  •  CH 

Ex  =  Ex  +Ex  +  Ex  +Ex 


(3) 


Although  exergy  is  extensive  property,  it  is  often  convenient  to  work  with  it  on  a  unit  of 
mass  or  molar  basis.  The  total  specific  exergy  on  a  mass  basis  may  be  written  as  follows: 


ex  =  exPH  +  exKN  +  exPT  +  exCH, 


(4) 


The  general  exergy  balance  can  be  written  as  follows: 

y  ^  Exin  y  ^  Exout  —  y  ^  Ex^t 


(5a) 


or 


Exheat  Ex 


work 


+  Ex 


mass,m 


—  Ex 


mass, out 


=  Exd 


est 


(5b) 


with 


Exheat  =  211- 


T 


0 


Ti 


Q 


k’ 


(6a) 


Exwop^  —  W , 


Fx 

J—/A. 


massjn 


=  EmoutiA 


out’ 


Exmass  out  —  X/h0ut*A 


out’ 


(6b) 

(6c) 

(6d) 


where  Qk  is  the  heat  transfer  rate  through  the  boundary  at  temperature  Tk  at  location  k 
and  W  is  the  work  rate. 

The  flow  (specific)  exergy  is  calculated  as  follows: 

i // =  (h  -  ho)  -  To(s  -  s0),  (7) 

where  h  is  enthalpy,  s  is  entropy,  and  the  subscript  zero  indicates  properties  at  the 
restricted  dead  state  of  P0  and  T0. 

The  rate  form  of  the  entropy  balance  can  be  expressed  as 
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where  the  rates  of  entropy  transfer  by  heat  transferred  at  a  rate  of  Qk  and  mass  flowing  at  a 
rate  of  m  are  Sheat  =  Qk/^k  and  Sm ass  =  ms,  respectively. 

Taking  the  positive  direction  of  heat  transfer  to  be  to  the  system,  the  rate  form  of  the 
general  entropy  relation  given  in  Eq.  (8)  can  be  rearranged  to  give 


^gen  —  ^^ouAout  ^^in^in  ^ 


Q 


k 


Ti 


Also,  it  is  usually  more  convenient  to  find  Sgen  first  and  then  to  evaluate  the  exergy 
destroyed  or  the  irreversibility  rate  /  directly  from  the  following  equation,  which  is  called 
Gouy-Stodola  relation  [22]: 


^O^gen- 


The  specific  exergy  (flow  exergy)  of  an  incompressible  substance  (i.e.,  water)  is  given 
by  [22] 


T  -  T0  -  T0  In 


The  total  flow  exergy  of  air  is  calculated  from  [23] 


<Aa,t  =  (Cp,a  +  coCp,w)T0[(T/T o)  -  1  -  In (T/T0)\  +  (1  +  1.6078co)i?ar0  In (P/P0) 


+  ^a^O 


(1  +  1.6078a;)  ln[(l  +  1.6078co0)  \ 
(1  +  1.6078a;)]  +  1.6078a;  \n(co/cD0)  / 5 


where  the  specific  humidity  ratio  is 


Assuming  air  to  be  a  perfect  gas,  the  specific  physical  exergy  of  air  is  calculated  by  the 
following  relation  [24] 


a, per 


\T  -  To  -  T0  In 


+  o  In 


p 

V 


3.2.  Energy  and  exergy  efficiencies 

Numerous  ways  of  formulating  exergetic  (or  exergy  or  second-law)  efficiency 
(effectiveness,  or  rational  efficiency)  for  various  energy  systems  are  given  in  detail 
elsewhere  [20].  It  is  very  useful  to  define  efficiencies  based  on  exergy  (sometimes  called 
Second  Law  efficiencies).  Whereas  there  is  no  standard  set  of  definitions  in  the  literature, 
two  different  approaches  are  generally  used — one  is  called  1 " ‘brute-force”,  while  the  other 
is  called  “functional.”  [25]. 

•  A  “brute-force”  exergy  efficiency  for  any  system  is  defined  as  the  ratio  of  the  sum  of  all 
output  exergy  terms  to  the  sum  of  all  input  exergy  terms. 

•  A  “functional”  exergy  efficiency  for  any  system  is  defined  as  the  ratio  of  the  exergy 
associated  with  the  desired  energy  output  to  the  exergy  associated  with  the  energy 
expended  to  achieve  the  desired  output. 
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Here,  in  a  similar  way,  exergy  efficiency  is  defined  as  the  ratio  of  total  exergy  output  to 
total  exergy  input,  i.e. 


(15) 


It  is  clear  that  the  brute-force  definition  can  be  applied  in  a  straightforward  manner, 
irrespective  of  the  nature  of  the  component,  once  all  exergy  flows  have  been  determined. 
The  functional  definition,  however,  requires  judgment  and  a  clear  understanding  of  the 
purpose  of  the  system  under  consideration  before  the  working  equation  for  the  efficiency 
can  be  formulated  [25]. 

3.3.  Exergetic  improvement  potential 

Van  Gool  [26]  has  also  proposed  that  maximum  improvement  in  the  exergy 
efficiency  for  a  process  or  system  is  obviously  achieved  when  the  exergy  loss  or 
irreversibility  (Eiin  —  Exout)  is  minimized.  Consequently,  he  suggested  that  it  is  useful 
to  employ  the  concept  of  an  exergetic  “improvement  potential”  when  analyzing  different 
processes  or  sectors  of  the  economy.  This  improvement  potential  in  the  rate  form,  denoted 
IP,  is  given  by 


IP  =  (1  -  e)(Exm  -  Exout). 


(16) 


3.4.  Some  thermodynamic  parameters 

Thermodynamics  analysis  of  RERs  may  also  be  performed  using  the  following 
parameters  [27]: 

Fuel  depletion  ratio: 


Relative  irreversibility: 


(18) 


Productivity  lack: 


Exergetic  factor: 


(19) 
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4.  Exergetic  analysis  and  evaluation  of  renewable  energy  resources 

4.1.  Solar  energy  systems 


In  evaluating  the  performance  of  solar  energy  systems  using  exergy  analysis  method, 
calculation  of  the  exergy  of  radiation  is  very  crucial.  However,  its  calculation  is  a  problem 
of  unquestionable  interest,  since  exergy  represents  the  maximum  quantity  of  work  that  can 
be  produced  in  some  given  environment  (usually  the  terrestrial  environment,  considered  as 
an  infinite  heat  source  or  sink)  [28].  Over  a  period  of  more  than  20  years,  many  papers 
including  various  approaches  to  this  calculation  have  been  published,  as  reported  in 
[29,30].  Among  these,  the  first  one  [31],  which  is  called  by  Petela  [29,32]  the  maximum 
efficiency  ratio  (or  exergy-to-energy  ratio  for  radiation)  for  determining  an  exergy  of 
thermal  emission  at  temperature  T,  is  as  follows: 


,max 


4  To 
3  T  ’ 


where  T  was  taken  to  equal  the  solar  radiation  temperature  (Ts)  with  6000  K  in  exergetic 
evaluation  of  a  solar  cylindrical-parabolic  cooker  and  a  solar  parabolic-cooker  by  Petela 
[32]  and  Ozturk  [33],  respectively.  Eq.  (21)  was  also  derived  by  Szargut  [34]  who  presented 
a  simple  scheme  of  a  reversible  prime  mover  (different  from  that  of  Petela)  transforming 
the  energy  of  radiation  into  mechanical  or  electrical  work  and  took  into  account  that  the 
solar  radiation  had  a  composition  similar  to  that  of  a  black  body. 

The  history  of  Eq.  (21)  is  as  follows  [35]:  It  was  first  derived  by  Petela  [31]  and 
independently  by  Landsberg  and  Mallinson  [36]  and  Press  [37],  while  Spanner  [38] 
proposed  Petela’s  equation  for  the  direct  solar  radiation  exergy  [29].  As  a  consequence,  it  is 
sometimes  called  the  Landsberg  efficiency  [34],  the  Petela-Landsberg  upper  bound  [39],  or, 
more  appropriately,  the  Petela-Press-Landsberg  and  Mallinson  formula  [28]. 

Petela  [29]  reported  that  in  the  existing  literature  in  the  above  field,  the  most  discussed 
formulae  for  heat  radiation  exergy  were  those  derived  by  him  [31],  Spanner  [38]  and  Jeter 
[40]  given  in  Eqs.  (21),  (24b),  and  (24c),  respectively.  As  for  as  the  application  of  these 
equations  to  various  thermal  systems  by  some  investigators,  Chaturvedi  et  al.  [41]  used 
Eq.  (21)  in  the  performance  assessment  of  a  variable  capacity  direct  expansion  solar- 
assisted  heat  pump  system  for  domestic  hot  water  application.  Eskin  [42]  considered 
Eq.  (24b)  in  evaluating  the  performance  of  a  solar  process  heating  system,  while  Singh  et 
al.  [12],  and  Ucar  and  Inalli  [43]  considered  Eq.  (24c)  in  the  exergetic  analysis  of  a  solar 
thermal  power  plant  and  the  exergoeconomic  analysis  of  a  solar-assisted  heating  system, 
respectively. 

Petela  [29]  also  listed  three  various  limiting  energy  efficiencies  of  utilization  of  the 
radiation  matter  in  Table  1  [29]  and  explained  them  as  follows:  “Bejan’s  conclusion  was 
that  all  theories  (Spanner’s,  Jeter’s  and  Petela’s)  concerning  the  ideal  conversion  of  thermal 
radiation  into  work,  although  obtaining  different  results,  were  correct.  Although  really 
correct,  they  had  incomparable  significance  for  the  true  estimation  of  the  radiation  exergy 
value.  All  of  the  efficiencies  assumed  work  as  an  output.  However,  Petela’s  work  was  equal 
to  the  radiation  exergy,  Jeter’s  work  was  the  heat  engine  cycle  work  and  Spanner’s  work 
was  an  absolute  work.  As  an  input,  Petela  and  Spanner  assumed  the  radiation  energy, 
whereas  Jeter  assumed  heat  [29].” 
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Table  1 

Numerators  (output)  and  denominators  (input)  of  the  limiting  energy  efficiency  of  radiation  utilization  by  three 
different  researchers  [29] 


Investigators 

Input 

Output 

Unified  efficiency  expression 

Petela  [31] 

Radiation  energy 

Useful  work  =  radiation  exergy 

i  +  ip°Y-4r° 

3 \T  J  3  T 

(210 

Spanner  [38] 

Radiation  energy 

Absolute  work 

x  4  T0 

3  T 

(24b) 

Jeter  [40] 

Heat 

Net  work  of  a  heat  engine 

i-r» 

T 

(24c) 

Svirezhev  et  al.  [44]  applied  a  new  concept  regarding  the  exergy  of  solar  radiation  to  the 
analysis  of  satellite  data  describing  the  seasonal  (monthly)  dynamics  of  four  components  of 
the  global  radiation  balance:  incoming  and  outgoing  short  and  long-wave  solar  radiation. 
They  constructed  maps  of  annual  mean  exergy  and  some  derivative  values  using  real 
satellite  data.  They  concluded  that  exergy  could  be  a  good  indicator  for  crucial  regions  of 
the  oceans,  which  were  indicated  where  its  values  were  maximal. 


4.1.1.  Solar  collector  applications 

The  instantaneous  exergy  efficiency  of  the  solar  collector  can  be  defined  as  the  ratio  of 
the  increased  water  exergy  to  the  exergy  of  the  solar  radiation  [33].  In  other  words,  it  is  a 
ratio  of  the  useful  exergy  delivered  to  the  exergy  absorbed  by  the  solar  collector  [12]. 

_  Exu 

^scol  —  v:  V^^-7 


Ex 


scol 


with 


or 


Exu  —  fi^w[(^w,out  ^wfm)  T()(£W5out  ^w,in)] 


Exu  =  mwC 


w 


(T’w.out  -  T’w.in)  -T(A  Ln 


T 


w,out 


T 


w,m 


Exu  =  Q 


U 


1  - 


T 


T  w,out  ^w,in 


Lll 


T 


w,out 


T 


w,m 


and 


EXscol  —  zl/plAsrad, 


max’ 


(23a) 


(23b) 


(23c) 


(24a) 


where  the  subscripts  “in”  and  “out”  denote  the  inlet  and  outlet  of  the  solar  collector, 
respectively,  while  *Asrad  max  maY  a^o  be  calculated  using  Eqs.  (21),  (24b)  and  (24c)  given  in 
Table  1  [29]. 

The  exergy  of  solar  radiation  with  beam  (/be)  and  diffuse  (/d)  components  for  parabolic 
collectors  is  also  given  by  the  following  equation  [45],  which  was  used  by  Eskin  [46]  in  the 
performance  evaluation  of  a  cylindrical  parabolic  concentrating  collector. 


*Asrad,be,d  ^be 


4  T 


o 


372 


+  ^d (  1  — 


4T 


o 


3  T 


* 


(25) 
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with 

jS  =  0.9562  +  0.2777  ln(l /k)  +  0.051  Ik,  (26) 

^  S 

where  k  is  the  dilution  factor  of  diffuse  radiation  [47]. 


4. 1.1.1.  Solar  water  heating  systems.  Solar  water  heater  is  the  most  popular  means  of 
solar  energy  utilization  because  of  technological  feasibility  and  economic  attraction 
compared  with  other  kinds  of  solar  energy  utilization.  Fig.  1  illustrates  a  schematic  of  a 
typical  domestic-scale  solar  water  heater,  which  mainly  consists  of  a  solar  collector  and  a 
storage  barrel  [48]. 

Exergy  efficiency  of  the  domestic-scale  solar  water  heater  may  be  calculated  using 
Eq.  (15)  as  follows: 


Ex 


^s, heater  — 


output 


Ex 


(27) 


J-^sun 

where  exergy  from  the  storage  barrel  to  the  end-user  (output  exergy)  may  be  written  as 
follows  [48],  assuming  the  temperature  distribution  in  the  storage  barrel  is  linear: 


Ex 


output 


—  rhxtj  C 


Ttoy  +  T 


bottom 


w  ^  w 


^bottom  T^m  w  Cw 


In 


mw  Cw  T o  (  In 


T 


top 


T 


-  1 


o 


T  top  T  bottom  T  bottom 


(28) 


and  Exsun  is  the  exergy  from  sun  (input  exergy)  and  may  be  found  using  Eq.  (24a). 

The  collector  exergy  efficiency  is  calculated  by  the  following  equation,  while  it  may  be 
found  using  Eqs.  (22),  (23b)  and  (24). 


Ex 


£s,col  — 


col, tank 


Ex 


(29) 


sun 


where  Excoi,tank  is  exergy  from  the  collector  to  the  storage  barrel  and  may  be  written  in  a 
similar  way  given  by  Eq.  (23b)  as  follows  [48]. 


Ex 


col, tank 


=  mw  C 


w  ^  w 


( -Ew  out  —  To)  —  To  ( In 


T 


w,out 


T 


0 


(30) 


Fig.  1.  A  typical  domestic-scale  solar  water  heater  [48]. 
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Xiaowu  and  Ben  [48]  performed  an  exergy  analysis  of  a  domestic-scale  solar  heater.  They 
used  the  data  obtained  from  Domestic-Scale  Solar  Water  Heaters  Performance  Testing 
Station,  located  in  Yun  Nan  province,  China  as  follows:  Top  and  bottom  hot  water 
temperatures  of  50  and  38  °C  in  the  storage  barrel,  respectively,  2.5  m  of  collector  area,  an 
average  ambient  temperature  of  25  °C,  an  average  ambient  water  temperature  of  24  °C, 
466  W/m  of  whole-day  average  solar  radiation  and  196.4  kg  of  water  storing  log  capacity 
for  the  storage  barrel.  They  calculated  the  exergy  efficiency  value  of  the  solar  water  heater 
using  the  above  data  and  found  to  be  0.77%.  The  main  conclusions  drawn  from  the  results 
of  their  study  were  as  follows:  The  exergy  efficiency  of  domestic-scale  water  heater  was  small 
as  the  output  energy  was  of  low  quality.  The  storage  barrel  had  large  exergy  losses.  It  was 
worthwhile  to  design  a  new  style  of  storage  barrel  for  reducing  unnecessary  mixing  of  water. 
In  order  to  improve  the  exergy  efficiency  of  domestic-scale  water  heater,  a  judicious  choice 
of  width  of  plate  and  layer  number  of  cover  was  necessary  [48]. 

4. 1.1.2.  Solar  space  heating  and  cooling 

4. 1.1. 2.1.  Air  systems.  Air  systems  are  indirect  water  heating  systems  that  circulate  air 
via  ductwork  through  the  collectors  to  an  air-to-liquid  heat  exchanger.  In  the  heat 
exchanger,  heat  is  transferred  to  the  potable  water,  which  is  also  circulated  through  the 
heat  exchanger  and  returned  to  the  storage  tank.  A  double  storage  tank  system  is  shown  in 
Fig.  2  [48].  This  type  of  system  is  used  most  often,  because  air  systems  are  generally  used 
for  preheating  domestic  hot  water  and  thus  auxiliary  is  used  only  in  one  tank.  The  main 
advantage  of  the  system  is  that  air  does  not  need  to  be  protected  from  freezing  or  boiling,  is 
noncorrosive,  and  is  free.  The  disadvantages  are  that  air  handling  equipment  (ducts  and 
fans)  need  more  space  than  piping  and  pumps,  air  leaks  are  difficult  to  detect,  and  parasitic 
power  consumption  is  generally  higher  than  that  of  liquid  systems  [49]. 
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Fig.  2.  A  schematic  of  air  system  with  a  double  storage  tank  drawn  by  Xiaowu  and  Ben  [48]. 
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The  relations  used  to  evaluate  air  systems  using  exergy  analysis  method  are  presented  by 
the  author  since  any  exergetic  assessment  based  on  the  actual  data  could  not  be  found.  The 
main  parts  to  be  included  in  the  exergy  analysis  are  (I)  air  solar  collector,  (II)  air-to-water 
hat  exchanger,  (III)  fan,  (IV)  storage  tank  and  (V)  circulating  pump.  The  exergy 
destruction  and  exergy  efficiency  relations  for  each  component  of  the  whole  air  system  are 
given  below  using  Eqs.  (5a)  and  (15),  respectively. 

I.  Air  solar  collector : 


Exdestj  —  Exjn?air  Ex 


out, air 


+  Ex 


s? 


A 


Exout,air  Ex 


m,air 


Exscoi 


where  Exin?air  and  Exout,air  are  calculated  using  Eqs.  (6c),  (6d),  (12)  and  (14),  while 
Exscoi  is  found  from  Eq.  (24a). 

II.  Air-to-heat  exchanger : 


EXdest,II  —  (EXh0t  , air, in  T  Excoid,w,in  )  -  (Exhot  , air, out  T  Excoid,w,out); 


fill 


Ex 


cold, w, out 
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III.  Fan\ 
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(36) 


(37) 


IV.  Storage  tank : 

Exdestjv  and  si\  may  be  calculated  in  a  similar  way  given  by  Eqs.  (33)  and  (34). 

V.  Circulating  pump : 

EXdest,VI  =  ^AvOAin  —  4* out, a)  "E  fEpump,e?  (38) 


8\  = 


Exout9a  EXin  ^hw(*Aout,a  *Ain) 
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pump,e 


with 


fEpump,e  lEpump  /  id  pump,e  h  pump,mech  )  * 
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where  t]pum p,e  and  ^pump,mech  are  the  circulating  pump  motor  electrical  and  circulating 
pump  mechanical  efficiencies. 


The  overall  exergy  destruction  in  the  whole  air  system  is  calculated  as  follows,  while  the 
overall  exergy  efficiency  is  found  on  a  product/fuel  basis  given  by  Eq.  (15). 

EXdest5o  —  EXdestj  +  EXdestji  T  EXdestjn  T  EXdestjv  T  Ex^es^y.  (41) 


4.1. 1.2.2.  Water  systems.  There  are  many  variations  of  systems  used  for  solar  space 
heating.  The  exergetic  analysis  of  these  systems  may  be  performed  using  relations  given  in 
Sections  4. 1.1.1  and  4. 1.1. 2.1. 

4. 1.1. 2. 3.  Heat  pumps.  Heat  pumps  are  preferred  and  widely  used  in  many 
applications  due  to  their  high  utilization  efficiencies  compared  to  conventional  heating 
and  cooling  systems.  There  are  two  common  types  of  heat  pumps:  air-source  heat  pumps 
and  ground-source  (geothermal)  heat  pumps  (GSHPs).  The  main  components  of  these 
systems  are  a  compressor,  a  condenser,  an  expansion  valve  and  an  evaporator.  In  addition, 
circulating  pumps  and  fans  are  used.  GSHPs  are  exergetically  studied  in  more  detail  in 
Section  4. 3. 2. 2  and  relevant  exergy  relations  are  presented  in  a  tabulated  form,  while 
studies  conducted  by  various  investigators  are  compared.  The  relations  given  in  these 
sections  may  be  used  in  modeling  and  exergetic  evaluation  for  heat  pump  systems.  In  case 
that  a  solar  assisted  heat  pump  with  direct  expansion  of  the  refrigerant  within  the  solar 
collector  is  used,  the  collector-evaporator  part  may  be  analyzed  using  Eq.  (24a)  for  solar 
radiation  input. 

Tsaros  et  al.  [50]  studied  on  exergy  analysis  of  heat  pumps.  Exergy  transports  between 
the  components,  to  the  loads  and  to  the  environment  were  assessed.  The  exergy 
consumption  in  each  component  were  also  determined.  Typical  results  were  presented  for 
the  operation  of  a  nominal  3.5  ton  residential  air-to-air  heat  pump,  in  both  the  heating 
and  cooling  modes.  Exergy  efficiency  values  for  various  system  components  are  listed  in 
Table  2  [50].  Of  the  major  components,  the  outdoor  coil  had  the  highest  exergy  efficiency, 
followed  by  the  compressor  and  the  indoor  coil.  The  exergy  efficiency  of  the  indoor  coil 
varied  greatly  with  load  and  exceeded  that  of  the  compressor  under  some  conditions.  The 
fans,  with  the  values  less  than  5%,  were  the  least  efficient  [50]. 

Bilgen  and  Takahashi  [51]  carried  out  exergy  analysis  of  heat  pump — air  conditioner 
systems.  The  experimental  system,  shown  in  Fig.  3  [51],  was  Matsushita  room  air 
conditioner  CS-XG28M.  The  irreversibilities  due  to  heat  transfer  and  friction  were 
considered.  The  coefficient  of  performance  based  on  the  first  law  of  thermodynamics  as  a 
function  of  various  parameters,  their  optimum  values,  and  the  efficiency  and  coefficient  of 
performance  based  on  exergy  analysis  were  derived.  Based  on  the  exergy  analysis,  a 
simulation  program  was  developed  to  simulate  and  evaluate  experimental  systems.  It  was 
found  that  the  exergy  efficiency  varied  from  37%  to  25%  both  a  decreasing  function  of 
heating  or  cooling  load.  The  exergy  destructions  in  various  components  were  also 
determined  for  further  study  and  improvement  of  its  performance.  It  was  suggested  that  to 
improve  the  performance  of  the  heat  pump  system,  each  component  could  be  further 
studied  from  exergy  usage  and  economics  points  of  view.  Eventually,  a  system  exergy- 
economics  optimization  study  could  be  carried  out  to  identify  each  component’s  relative 
importance  with  respect  to  operational  conditions. 
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Table  2 

Exergy  efficiency  values  of  the  system  components  [50] 


Component 

Exergy  efficiencies  in 

%  at  various  outdoor  temperatures 

-18.75  °C  (17°F) 

8.33  °C  (47°F) 

35  °C  (95°F) 

Indoor  coil 

82.6 

57.0 

35.3 

Compressor,  4-way  valve,  outdoor  coil,  throttle 

36.0 

40.8 

29.8 

Compressor,  4-way  valve,  outdoor  coil 

44.0 

47.7 

35.2 

Throttle 

81.6 

85.5 

84.5 

Compressor,  4-way  valve 

47.9 

55.0 

53.7 

Outdoor  coil 

92.0 

86.9 

65.6 

Compressor 

50.5 

58.5 

58.3 

4-way  valve 

94.7 

93.9 

92.1 

Outdoor  fan  and  motor 

4.41 

4.41 

4.41 

Indoor  fan  and  motor 

3.49 

3.49 

5.81 

Overall  without  fans 

30.1 

23.4 

10.7 

Overall 

21.7 

18.1 

8.43 

Compressor  includes  associated  hardware,  while  4-way  valve  includes  piping  to  coils. 


EXPANSION  VALVE 


Fig.  3.  Schematic  of  the  experimental  set  up  showing  measurement  system  and  instrumentation  drawn  by  Bilgen 
and  Takahashi  [51]. 

Aprea  et  al.  [52]  studied  on  the  performances  of  a  vapor  compression  experimental 
system  working  both  as  a  water  chiller  and  as  a  heat  pump,  using  as  refrigerant  fluids  R22 
and  its  substitute  R417A  (R125/R134a/R600  46.6/50/3.4%  in  mass).  This  kind  of  system 
has  been  used  in  industrial  processing  or  air  conditioning  plants  where  a  supply  of 
refrigerated  and  heated  water  was  required.  The  experimental  system  was  principally  made 
up  of  a  unit  that  could  operate  both  as  refrigeration  system  and  heat  pump.  It  was  made  up 
of  a  hermetic  compressor  Scroll,  a  plate-type  water  heat  exchanger  inserted  in  a  water  tank, 
a  finned  tube  air  heat  exchanger,  two  thermostatic  expansion  valves. 
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Aprea  et  al.  [52]  evaluated  the  efficiency  defect  for  each  component  of  the  system, 
considering  the  ratio  between  the  exergy  destroyed  in  each  one  and  the  exergy  required  to 
sustain  the  process  (i.e.,  the  electrical  power  supplied  to  the  compressor).  The  efficiency 
defects  of  R417A  were  always  higher  than  the  defects  of  R22  of  about  6%  related  to  the 
valve  and  the  heat  exchangers,  and  of  about  10%  referring  to  the  compressor.  Moreover,  it 
was  experimentally  verified  that  the  COP  of  R22  was  higher  than  that  of  R417A  also  of 
about  18%  related  to  the  water  chiller  system  and  of  about  15%  referring  to  the  heat 
pump. 

Ma  and  Li  [53]  exergetically  analyzed  a  heat  pump  system  with  economizer  coupled  with 
scroll  compressor  illustrated  in  Fig.  4  [53]  and  derived  expressions  for  exergy  losses  of  the 
processes  in  the  system,  exergy  loss  ratio  and  efficiency  of  the  system.  They  also  determined 
the  exergetical  losses  and  efficiency  of  this  heat  pump  system  using  experimental  data, 
while  they  assessed  characteristics  of  exergy  variation  and  transport  and  exergetic 
efficiency  of  the  system.  The  experimental  system  studied  is  explained  in  more  detail 
elsewhere  [54]. 

Using  Eq.  (15),  Ma  and  Li  [53]  expressed  exergy  efficiency  as  a  ratio  of  the  effective 
exergy  output  of  the  heat  pump  (equal  to  the  heat  load  exergy  in  the  condenser  plus  the 
cool  capacity  exergy  in  the  evaporator)  to  the  exergy  of  the  electrical  energy  as  follows: 

_  E^eff, output  _  Exh  +  EXcool 

£  Fy  Ex 

From  the  distribution  of  the  exergy  loss,  it  was  found  that  the  maximum  loss  occurred  in 
the  compressor,  accounting  for  about  82%  of  the  total  exergy  destruction.  This  value 
would  increase  when  the  temperature  difference  between  the  condensation  and  the 
evaporation  could  be  increased.  The  compressor  was  a  bottleneck  for  the  efficiency  of  the 
heat  pump  system,  and  increasing  energy  efficiency  ratio  for  the  compressor  was  crucial  to 
the  improvement  of  the  heat  pump  system.  For  hermetic  compressor,  there  were  two 
mechanisms  for  exergy  to  consumption:  one  is  was  the  process  in  which  electrical  power 
was  converted  to  mechanical  work  in  motor,  and  the  other  was  the  process  of  the 
mechanical  work  being  converted  to  the  pressure  energy  of  the  refrigerant  in  the 
compressor.  The  exergy  losses  in  the  evaporator  and  condenser  were  found  to  be  8.9%  and 


scroll  compressor 


condenser 


Fig.  4.  Schematic  of  a  heat  pump  system  with  economizer  coupled  with  scroll  compressor  drawn  by  Ma 
and  Li  [53]. 
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5.3%  of  the  total  loss,  respectively,  and  were  10.8%  and  6.5%  of  the  loss  in  the 
compressor.  Although  they  were  relatively  small  compared  with  the  loss  in  the  compressor, 
minimizing  the  exergy  losses  in  evaporator  and  condenser  was  also  important  for  the 
improvement  of  the  efficiency  of  the  heat  pump,  and  improving  the  evaporator  was  more 
beneficial  for  the  overall  system  efficiency  than  improving  the  condenser.  The  exergy  losses 
existing  in  the  essential  components  of  the  heat  pump  system  such  as  the  compressor,  the 
evaporator  and  the  condenser  were  predominant  in  the  total  exergy  loss,  namely  96.2%  of 
the  total.  In  summary,  the  best  gain  in  the  heat  pump  system  performance  could  be 
obtained  through  improvement  in  the  compressor,  followed  by  the  evaporator  and 
condenser.  Using  Eq.  (42),  exergy  efficiency  of  the  air-source  heat  pump  system  was 
determined  to  vary  approximately  from  33%  to  42%  at  the  evaporation  temperatures  of 
—25  to  — 12  °C  with  a  fixed  condensing  temperature  of  45  °C  [53]. 

Badescu  [55]  proposed  a  model  for  the  heating  system  of  an  ecological  building 
whose  main  energy  source  was  solar  radiation.  The  most  important  component  of  the 
heating  system  was  a  vapor  compression  heat  pump.  The  building’s  electric  energy  was 
provided  by  a  photovoltaic  (PV)  array  or  from  the  AC  grid,  when  necessary.  Exergy 
analysis  of  the  heat  pump  was  performed.  The  state  parameters  and  the  process  quantities 
were  evaluated  by  using,  as  input,  the  building  thermal  load.  It  was  reported  that  most  of 
the  exergy  losses  occurred  during  compression  and  condensation.  Preliminary  results 
indicated  that  the  PV  array  could  provide  all  the  energy  required  to  drive  the  heat  pump 
compressor,  if  an  appropriate  electrical  energy  storage  system  was  provided.  In  the 
analysis,  meteorological  data  measured  by  the  Romanian  Meteorological  and  Hydro- 
logical  Institute  during  1961  in  Bucharest  (latitude  44.5°N,  longitude  26.2°E,  altitude 
131m  above  sea  level)  were  used,  while  the  exergy  efficiency  of  the  heat  pump,  was  found 
using  the  following  equation: 


£hp 


Wcomp,C 

Wcomp,a 


where  wcomp,a  and  wC0mp,c  are  the  specific  mechanical  works  entering  the  actual  and 
reversible  (Carnot)  compressor,  respectively.  Using  Eq.  (43),  exergy  efficiency  of  the  heat 
pump  was  found  to  be  15.53%  and  18.37%  for  the  refrigerants  of  R114  and  R12, 
respectively,  based  on  the  results  obtained  in  the  cold  day  of  January  22,  1961  at  6.00  a.m. 
in  Bucharest,  when  the  ambient  temperature  was  255.7  K  [55]. 

Cervantes  and  Torres-Reyes  [56]  performed  an  experimental  study  of  a  solar  assisted 
heat  pump  with  direct  expansion  of  the  refrigerant  within  the  solar  collector  and  evaluated 
the  performance  of  this  system  using  exergy  analysis  method.  The  condensation  of  the 
refrigerant  took  place  in  an  R22/air  condenser  made  out  of  stainless-steel,  with  a  heat 
transfer  area  of  9.8  m2.  The  compression  stage  was  carried  out  by  a  hermetic  compressor, 
and  the  expansion  of  R-22  occurred  in  a  thermal  expansion  valve.  The  refrigerant  was 
evaporated  inside  an  uncovered  solar  flat  collector,  with  a  collection  area  of  364.5  m  . 
Experimental  tests  were  performed  at  a  solar  radiation  range  of  200-1 100  W/m  and 
ambient  temperatures  varying  from  20  to  32  °C.  The  actual  amount  of  work  supplied  to 
the  compressor,  given  by  the  electricity  consumption,  ranged  from  1.1  to  1.36  kW.  The 
maximum  exergy  efficiency,  defined  as  the  ratio  of  the  outlet  to  the  inlet  exergy  flow 
in  every  component  of  the  heat  pump  cycle,  was  determined  taking  into  account  the 
typical  parameters  and  performance  coefficients.  The  results  of  this  exergy  analysis 
pointed  out  that  the  main  source  of  irreversibility  could  be  found  in  the  evaporator  of 
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the  heat  pump  (i.e.,the  solar  collector)  emphasizing  that  incoming  solar  radiation,  as 
shown  in  Fig.  5  [56].  A  full  analysis  of  the  inlet  and  outlet  exergy  flows  and  the  main 
parameters  associated  to  them,  for  each  component  of  the  solar  assisted  heat  pump,  is 
given  in  Ref.  [57]. 

Ucar  and  Inalli  [43]  developed  an  exergoeconomic  model  for  analysis  and  optimization 
of  solar-heating  systems  with  residential  buildings.  The  model  system  consisted  of  flat  plate 
solar  collectors,  a  heat  pump,  an  underground  storage  tank  and  a  heating  load,  as  shown 
in  Fig.  6  [43].  The  heat  produced  by  the  solar  collectors  was  stored  in  underground  storage 
tanks  throughout  the  year.  This  heat  was  only  extracted  during  the  heating  season.  It  was 
found  from  the  exergetic  analysis  results  for  the  cylindrical  tank  system  that  the  largest 
exergy  losses  occurred  in  the  house  at  78.7%,  followed  by  the  heat  pump  at  40.9%,  the 
collector  at  39.7%  and  the  storage  at  19.8%. 


Fig.  5.  Grassmann  (exergy  flow)  diagram  of  the  solar  assisted  heat  pump  system  with  direct  expansion  of  the 
refrigerant  within  the  solar  collector  drawn  by  Cervantes  et  al.  [56]. 


SOLAR  COLLECTOR 


Fig.  6.  Schematic  diagram  of  a  solar-assisted  heating  system  for  residential  buildings  drawn  by  Ucar  and  Inalli 
[43], 
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4.1.13.  Solar  refrigeration 

4.1.13.1.  Solar-driven  ejector  refrigeration  systems.  A  solar-driven  ejector  refrigera¬ 
tion  system  consists  of  a  solar  collector  subsystem  and  an  ejector  refrigeration  subsystem. 
The  major  components  in  the  refrigeration  cycle  are  an  ejector,  a  condenser,  a  generator, 
an  evaporator,  an  expansion  device  and  a  pump.  Fig.  7  [58]  illustrates  a  schematic  of  this 
system,  which  was  energetically  analyzed  by  Pridasawas  and  Lundqvist  [58]  in  more  detail. 
In  assessing  this  system  using  exergy  analysis  method,  the  relations  given  in  Section  4. 3. 2. 2 
and  additionally  the  exergetic  efficiency  definitions  explained  below  may  be  used. 

The  exergy  inputs  and  outputs  may  take  different  forms,  e.g.,  radiation,  heat  or 
electricity.  Thus,  the  exergetic  efficiency  of  the  ejector  refrigeration  cycle  is  defined  as  the 
ratio  of  the  cooling  exergy  (at  the  evaporator)  and  the  exergy  inputs  to  the  generator  and 
the  pump: 


Ex, 


evap 


Exg  "E  Exp  e 


with 


Exg  =  eg(' (45) 

where  Q  is  energy  input  to  the  generator  and  Tg  is  the  generating  temperature  [58]. 

The  system  exergetic  efficiency  is  defined  as  the  ratio  of  the  cooling  exergy  and  the 
exergy  of  the  available  solar  heat  input  plus  the  electricity  input  to  the  pump: 


Ex 


evap 


Ex^h  +  Ekp5e 


with 


ava 


Solar  collector 


2  Generator 
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Condenser 


Throttle 
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Fig.  7.  Ejector  refrigeration  cycle  drawn  by  Pridasawas  and  Lundqvist  [58]. 
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where  gav a  is  available  heat  for  the  process  in  the  solar  collector  and  Tscol  is  the  average 
temperature  of  the  solar  collector  [58]. 

Pridasawas  and  Lundqvist  [58]  used  the  following  relation  established  by  De  Vos  [59]  for 
calculating  the  exergy  input  to  the  collector  per  unit  area  as  a  function  of  the  exergy 
emitted  from  the  sun  minus  the  albedo  of  the  earth  and  the  radiation  emitted  from  the 
solar  collector: 


Exscoi  =foT\  +  (1  -f)cT4plamt  -  CJT *co1,  (48) 

where  /  is  the  sunlight  dilution  factor,  which  equal  to  2.16  x  10-5  on  the  earth,  o  is 
Stefan-Boltzmann  constant  (5.67  x  10-8  W/m2  K4),  Ts  is  the  temperature  of  the  sun  (K) 
and  rplanet  is  the  temperature  of  the  planet  (K). 

The  exergy  loss  during  the  transformation  from  solar  radiation  to  heat  on  the  solar 
collector  is  given  as  [58]: 


/ 


scol,trans 


=  Exs  -  Exs?h. 


The  exergy  loss  from  the  exergy  input  to  the  solar  collector  to  the  working  fluid  can  be 
calculated  using  the  following  equations  [58]: 


with 


Exu  =  Qu  ( 1  -  Tl-V  (51) 

Pridasawas  and  Lundqvist  [58]  evaluated  the  performance  of  an  ejector  refrigeration 
cycle  driven  by  solar  energy  using  exergy  analysis  method.  The  analysis  was  based  on  the 
following  conditions:  a  solar  radiation  of  700  W/m  ,  an  evaporator  temperature  of  10  °C,  a 
cooling  capacity  of  5  kW,  butane  as  the  refrigerant  in  the  refrigeration  cycle  and  ambient 
temperature  of  30  °C  as  the  reference  temperature.  Irreversibilities  occurred  among 
components  and  depend  on  the  operating  temperatures.  The  most  significant  losses  in  the 
system  were  in  the  solar  collector  and  the  ejector.  The  latter  decreased  inversely 
proportional  to  the  evaporation  temperature  and  dominated  the  total  losses  within  the 
system.  For  the  above  operating  conditions,  the  optimum  generating  temperature  was 
about  80  °C.  The  exergy  losses  in  the  ejector  refrigeration  cycle  is  shown  in  Fig.  8  [58].  The 
highest  loss  with  about  50.5%  was  found  in  the  solar  collector,  followed  by  the  losses  in  the 
ejector,  which  was  15.8%  and  in  the  condenser  with  10.6%  (1.49 kW  =  0.89  kW  loss  and 
0.60  kW  rejected).  On  the  component  basis  of  the  ejector  refrigeration  subsystem,  the 
largest  loss  (32%)  occurred  in  the  ejector  followed  by  the  generator  (22%),  the  condenser 
(21%).  The  rest  of  the  exergy  loss  occurred  in  the  pump,  the  evaporator  and  the  expansion 
valve.  The  exergetic  efficiency  values  for  the  solar  collector  and  refrigeration  cycle  were 
14.53  and  4.39%,  respectively,  while  the  overall  exergetic  efficiency  was  0.66%. 

4. 1.1.3. 2.  Absorption  units.  Sozen  et  al.  [60]  investigated  a  prototype  of  an 
aqua-ammonia  absorption  heat  pump  system  using  solar  energy.  Fig.  9  [60]  illustrates  a 
schmatic  of  this  system.  The  aqua-ammonia  absorption  heat  pump  consisted  of  condenser, 
evaporator,  absorber,  generator,  pump,  expansion  valves,  refrigerant  heat  exchanger 
(precooler,  RHE)  and  mixture  heat  exchanger  (solution  heat  exchanger,  MHE).  When 
these  two  heat  exchangers  were  used  in  the  system,  the  performance  was  increased. 
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Solar  collector 
(4.25  kW) 
Condenser 
(1.49  kW) 
Ejector 
(1.33  kW) 
Generator 
(0.90  kW) 
Evaporator 
(0.18  kW) 

Pump 
(0.23  kW) 
Expansion  valve 
(0.04  kW) 


Fig.  8.  Exergy  balance  of  a  solar-driven  ejector  refrigeration  system  for  the  generating  temperature  of  90°C  and 
evaporation  temperature  of  10  °C  [58]. 


Fig.  9.  Schematic  of  a  absorption  heat  pump  drawn  by  Sozen  et  al.  [60]. 


For  each  individual  component  of  the  system  shown  in  Fig.  9,  the  exergy  destructions 
are  written  as  follows  [60]: 

Ex^estcond  =  Exg)  ^condt^  (-Eo/^cond)]? 

EXdest,evap  —  m7(Exio  -  Exn)  +  2ev apt1  -  (To/T  evap)]? 


(52) 

(53) 
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Eldest, absor  =  m7Exi2  +  m6Ex6  -  miEX]  -  gabsor[l  -  (To/Tateor)], 


EXdest,pump  — ril\(Ex.]  Ex2)  +  M'pump.e ? 


Exdest,g  =  m3Ex3  -  m4Ex4  -  to7Ex7  +  Qg[  1  -  ( T0/Tg )],  (56) 

Exdest,mix,HE  =  m2(Ex2  -  Ex3)  +  m4(Ex4  -  Ex5),  (57) 


Exdest,r,HE  =  w7(Ex8  -  Ex9  +  Exn  -  Exi2).  (58) 

The  exergetic  coefficient  of  performance  of  the  system  for  cooling  and  heating  are  as 
follows,  respectively  [60]: 

Scond  (  1 

copex,h  = - V 
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Sozen  et  al.  [60]  performed  the  performance  tests  of  the  system  under  the  climate 
condition  of  Ankara  in  Turkey.  A  parabolic  collector  was  installed  to  obtain  required 
temperatures.  In  the  experiments,  high  temperature  water  obtained  from  the  collector  was 
used  as  heat  source  needed  for  the  generator.  The  system  design  configuration  was 
analyzed  by  using  the  experimental  data.  The  effect  of  irreversibilities  in  thermal  process 
on  the  system  performance  in  the  system  were  determined.  Thermodynamic  analysis 
showed  that  both  losses  and  irreversibility  had  an  impact  on  absorption  system 
performance,  while  potential  for  the  improvement  in  the  components  of  the  system  was 
determined.  The  components  of  the  evaporator  and  absorber  of  the  system  had  a  higher 
exergy  loss  than  the  other  components.  In  order  to  improve  the  performance  of  the  system 
thermally,  those  components  should  be  modified.  Typical  COPex,h  values  of  the  system 
varied  from  0.13  to  0.45.  The  highest  and  lowest  COPex,h  values  were  obtained  at  higher 
evaporator  temperature  (10  °C)  and  the  generator  temperature  of  75  °C,  and  at 
revap  =  10  °C  and  Tcond  =  28  °C,  respectively. 

Izquierdo  et  al.  [61]  studied  on  the  entropy  generated,  the  exergy  destroyed  and  the 
exergetic  efficiency  of  lithium-bromide  absorption  thermal  compressors  of  single  and 
double  effect,  driven  by  the  heat  supplied  by  a  field  of  solar  thermal  collectors.  A  schematic 
of  a  single  effect  thermal  compressor  fed  by  a  thermal  collector  is  illustrated  in  Fig.  10  [61], 
where  the  solar  cases  heat  storage  system  was  not  included.  It  consisted  of  an  absorber,  a 
desorber,  the  heat  recoverer,  the  solution  pump  and  the  throttling  valve.  The  compressor, 
driven  by  the  heat  produced  by  a  field  of  solar  collectors,  absorbed  the  low  pressure. 

Izquierdo  et  al.  [61]  reported  as  follows:  Two  different  applications,  namely  air-cooled 
and  water-cooled  units,  were  taken  into  account  and  compared  with  each  other.  Water- 
cooled  compressors  operated  with  temperatures  and  pressures  lower  than  air-cooled 
compressors  considering,  in  both  cases,  the  same  suction  temperature,  equal  to  5  °C.  While 
the  absorption  temperature  in  water-cooled  compressors  could  reach  40  °C,  in  air-cooled 
systems  it  could  vary  between  30  °C  and  more  than  50  °C.  Under  these  conditions, 
the  discharge  temperature  (boiling  temperature  within  the  desorber)  of  a  single  effect 
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Fig.  10.  Schematic  of  a  single  effect  thermal  compressor  fed  by  a  thermal  collector  drawn  by  Izquierdo  et  al.  [61]. 


air-cooled  unit  ranged  from  65  to  110°C,  while  the  maximum  discharge  pressure  was 
around  0.1 2  bar.  The  discharge  temperatures  (boiling  temperature  within  the  high  pressure 
desorber)  of  the  double  effect  air-cooled  thermal  compressor  varied  from  1 10  °C  for  a  final 
absorption  temperature  of  30-180  °C  for  a  final  absorption  temperature  of  50  °C. 
Discharge  pressures  could  reach  values  of  0.3  and  1.5  bar,  respectively.  The  lithium- 
bromide  air-cooled  thermal  compressors  of  double-effect  could  be  viable  with  absorption 
temperatures  around  50  °C,  when  the  temperature  difference  between  the  lithium-bromide 
solution  and  the  outside  air  was  about  8  °C.  The  exergetic  efficiency  of  the  air-cooled 
compressor  of  double  effect  was  found  to  be  10%  greater  than  the  exergetic  efficiency  of 
the  air-cooled  thermal  compressor  of  single  effect.  The  exergetic  efficiency  of  the  air-cooled 
thermal  compressor  of  double  effect  was  also  determined  to  be  10%  higher  than  the 
exergetic  efficiency  of  the  water-cooled  thermal  compressor  of  single  effect.  The  exergy  of 
the  solar  radiation  was  better  used  by  the  air-cooled  thermal  compressor  as  the 
temperature  of  the  external  fluid  feeding  the  high  desorber  was  higher. 

4. 1.1.3. 3.  Adsorption  units.  Adsorption  refrigeration  is  considered  alternative  to  the 
conventional  vapor  compression  refrigerator,  especially  in  remote  areas  of  the  world 
without  grid  connected  electricity.  Activated  carbon/methanol,  activated  carbon/ammo¬ 
nia,  zeolite/water  and  silica  gel/water  are  the  adsorbent/adsorbate  pairs  commonly  used  in 
practically  realized  adsorption  refrigeration  cycles.  A  simple  adsorption  cycle  indicating  all 
the  heat  transfer  for  a  complete  cycle  is  shown  in  Fig.  11  [62].  Its  principle  of  operation 
involves  the  processes  of  isosteric  heating,  desorption,  isosteric  cooling  and  resorption  [62]. 

The  thermodynamic  efficiency  (exergetic  efficiency)  relates  the  actual  performance  of  a 
system  to  its  expected  maximum  performance  as  determined  from  its  ideal  Carnot  cycle.  It 
may  be  written  from  a  relation  proposed  by  Pons  et  al.  [63]  as  follows: 
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Fig.  11.  Schematic  of  a  simple  adsorption  refrigerator  [62], 


and 


COPa  =  ,  (63) 

2g> 

where  £>evap  and  are  the  heat  transferred  during  refrigeration  and  the  heat  used  to 
generate  refrigerant  during  generation,  respectively  [62]. 

Anyanwu  and  Ogueke  [62]  reported  as  follows:  The  COP  and  thermodynamic  efficiency 
depended  on  the  mass  of  adsorbate  generated  and  the  generation  temperature,  Tg.  Above  a 
certain  generator  temperature,  the  mass  of  adsorbate  generated  was  no  longer  a  direct 
function  of  the  temperature.  In  fact,  in  some  cases,  continued  increase  in  generator 
temperature  did  not  yield  any  adsorbate.  Thus,  an  improved  COP  and  consequently  the 
thermodynamic  efficiency  would  only  come  from  an  adsorbent  that  would  readily  desorb 
the  adsorbate  in  its  pores  up  to  a  low  level  of  about  20-25  g  of  adsorbate  per  kg  of 
adsorbent  for  the  range  of  temperature  obtainable  using  a  conventional  flat  plate  solar 
collector.  This,  however,  was  with  the  understanding  that  the  refrigerant  would  have  a 
high  latent  heat  of  vaporization.  The  activated  carbon/ammonia  for  0,-10  and  — 20  °C 
evaporator  temperatures  and  47  °C  condenser  temperature  gave  thermodynamic  efficiency 
values  above  100%  at  generator  temperatures  of  about  70  °C.  When  the  interest  was  air 
conditioning,  the  zeolite/water  pair  would  be  preferred  while  activated  carbon/ammonia 
would  be  the  preferred  pair  for  other  applications. 

4. 1.1.4.  Solar  cookers.  Solar  energy,  which  is  an  abundant,  clean  and  safe  source  of 
energy,  is  an  attractive  to  substitute  for  the  conventional  fuels  for  cooking  [33].  In  the 
literature,  there  are  two  studies  conducted  on  the  exergetic  evaluation  of  the  performance 
of  solar  cookers.  A  simply  designed  and  the  low-cost  parabolic-type  solar  cooker  (SPC) 
shown  in  Fig.  12  was  made  and  tested  by  Ozturk  [33].  He  experimentally  assessed  the 
energy  end  exergy  efficiencies  of  the  cooker.  Petela’s  analysis  [32]  was  inspired  by 
publication  of  Ozturk  [33],  in  which  for  the  first  time  the  exergy  efficiency  was  determined 
experimentally  for  the  SPC  of  the  cylindrical  trough  shape  [32]. 
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Fig.  12.  A  sketch  of  a  solar  parabolic  cooker  drawn  by  Ozturk  [33]. 


The  instantaneous  exergy  efficiency  of  the  solar  cooker  can  be  defined  as  the  ratio  of  the 
increased  water  exergy  to  the  exergy  of  the  solar  radiation.  Using  Eqs.  (21),  (22)  and  (23b), 
it  may  be  calculated  as  follows  [33]: 


m 


£scook,l  — 


w 


cw[(Tw,0ut  -  t\ v,in) -  r0(in(rWjOUt/rw>in))]/A? 

/TaCook[i  +  i/3(7yrs)4  -  (4/3)(r0/rs)] 


Exergy  efficiency  of  the  SPC  was  also  defined  by  Petela  [32]  as  the  ratio  of  the  exergy  of 
the  useful  heat  Q3u ,  at  temperature  T3,  and  of  the  exergy  of  solar  emission  at  temperature 
Ts.  Using  Eqs.  (21),  (22)  and  (24),  it  was  proposed  as  follows  [32]: 


£scook,2  — 


03u(i  -  (To/t3)) 


/t[i  +  l/3(T0/Ts)4  -  (4/3)(To/rs)] 
with  the  useful  heat  transferred  through  the  wall  of  the  cooking  pot  given  by 


(65) 


Q3u  =  -A3U3(T3-Tw),  (66) 

where  T3  is  the  absolute  temperature  of  the  surface  3,  A3  the  outer  surface  area  3  of  the 
cooking  pot,  U3  the  heat  transfer  coefficient  (which  takes  into  account  the  conductive  heat 
transfer  through  the  cooking  pot  wall  and  convective  heat  transfer  from  the  inner  cooking 
pot  surface  to  the  water)  and  Tw  the  absolute  temperature  of  water  in  the  cooking  pot,  as 
shown  in  Fig.  13  [32]. 

Fig.  12  illustrates  the  dimensions  of  a  SPC,  which  was  exergetically  tested  by  Ozturk 
[33].  The  overall  dimensions  of  the  SPC  are  0.90  x  1  m2.  The  SPC  was  constructed  of  steel 
profile  and  Cr-Ni  alloy  sheet.  Structural  frame  of  the  SPC,  made  of  steel  profile  of 
dimensions  0.15  x  0.15  m  and  Cr-Ni  alloy  sheet,  was  screwed  on  the  frame.  The  thickness 
of  Cr-Ni  sheet  was  0.0005  m.  In  the  center  of  the  concentrating  reflector,  a  cooking  pot  of 
0.10  m  width  and  0.05  m  in  depth  was  welded.  The  cooking  pot  was  made  of  galvanized 
steel  sheet  of  0.005  m  thickness.  All  outer  surfaces  of  the  cooking  pot  were  painted  with 
mat  black.  There  were  no  clouds  so  the  reflector  acted  only  on  direct  solar  radiation.  The 
reflector  could  be  adjusted  by  means  of  the  two  screws  on  both  sides  of  the  cooking  pot. 
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Fig.  13.  A  schematic  of  a  solar  parabolic  cooker  drawn  by  Petela  [32]. 


The  emissivity  of  the  cooking  pot  and  shield  surfaces  were  0.87  and  0.45,  respectively 
[33].  Assumed  were  the  following  data:  the  time  interval  A^  =  600s,  the  constant 
specific  heat  for  water  Cw  =  4.2kJ/kgK  and  the  solar  radiation  temperature  Ts  =  6000  K. 
The  experimental  time  period  was  from  10:00  to  14:00  solar  time.  During  this 
period,  it  was  found  that  the  daily  average  temperature  of  water  in  the  SPC  was  333  K 
and  the  daily  average  difference  between  the  temperature  of  water  in  the  cooking  pot 
and  the  ambient  air  temperature  was  31.6  K.  The  exergy  output  of  the  SPC  varied 
between  2.9  and  6.6  W.  The  exergy  efficiency  of  the  SPC  was  obtained  to  range  from  0.4% 
to  1.25%  [33]. 

Petela  [32]  derived  equations  for  heat  transfer  between  the  three  surfaces:  cooking  pot, 
reflector  and  imagined  surface  making  up  the  system.  His  model  allowed  for  theoretical 
estimation  of  the  energy  and  exergy  losses:  unabsorbed  insolation,  convective  and  radiative 
heat  transfer  to  the  ambient,  and  additionally,  for  the  exergy  losses:  the  radiative 
irreversibilities  on  the  surfaces,  and  the  irreversibility  of  the  useful  heat  transferred  to  the 
water.  The  exergy  efficiency  of  the  SPC  was  found  to  be  relatively  very  low  (about  1%), 
and  to  be  about  10  times  smaller  than  the  respective  energy  efficiency.  The  reasons  for  this 
were  mainly  due  to  the  large  exergy  of  the  escaping  insolation  and  additionally  due  to  the 
degradation  of  the  insolation  absorbed  on  the  surfaces  of  the  reflector  and  the  cooking  pot. 
The  theoretically  calculated  values  of  the  efficiencies  were  relatively  close  to  the  values 
experimentally  measured  by  Ozturk  [33].  The  influence  of  the  input  parameters 
(geometrical  configuration,  emissivities  of  the  surfaces,  heat  transfer  coefficients  and 
temperatures  of  water  and  ambience)  was  also  determined  on  the  output  parameters,  the 
distribution  of  the  energy  and  exergy  losses  and  the  respective  efficiencies  [32]. 

4. 1.1.5.  Industrial  process  heat.  Cylindrical  parabolic  collector  systems  look  very 
promising  for  delivering  industrial  process  heat  applications  in  95-350  °C  delivery 
temperature  range  [41].  In  this  regard,  Eskin  [42]  made  a  simulation  study  on  a  solar 
process  heat  system  shown  in  Fig.  14  [41],  validated  by  experimental  results  and 
represented  the  unsteady  performance  analysis  of  this  system.  She  also  evaluated  exergetic 


622 


A.  Hepbasli  /  Renewable  and  Sustainable  Energy  Reviews  12  ( 2008 )  593-661 


III 


Return 


Fig.  14.  A  schematic  view  of  the  solar  process  heating  system  [41]  (I:  cylindrical  parabolic  collector  array,  II: 
storage  tank,  III:  circulating  pump,  IV:  auxiliary  heater). 


performance  of  the  particular  systems,  while  her  analysis  considered  the  unsteady-state 
thermal  and  exergetic  assessment  of  the  collector  and  the  storage  tank  individually. 

The  exergy  efficiency  of  the  cylindrical  parabolic  solar  collector  was  calculated  as 
follows: 


£par,col  — 


|  7par,col 


Ex 


par,s 


with 


Ekpars  ^*Asrad,be,d*  (68) 

The  auxiliary  heater  in  the  system  eliminated  the  difference  between  energy  demand  and 
solar  energy  harvesting.  When  the  temperature  of  the  water  entering  the  load  fell  below  the 
predetermined  temperature  level,  the  auxiliary  heater  was  activated.  The  simulation  results 
indicated  that  the  exergetic  efficiency  (at  a  maximal  value  of  about  37%)  was  highly 
dependent  on  the  ratio  of  mass  flow  rates  and  the  use  of  an  auxiliary  heater  in  the  system. 


4. 1.1.6.  Solar  desalination  systems.  There  are  very  limited  studies  on  the  exergetic 
evaluation  of  solar  desalination  systems  in  the  open  literature.  The  Spanish  research 
institution,  Centro  de  Investigaciones  Energtticas,  Medioambientales  y  Tecnologicas  and 
the  German  Deutsche  Forschungsanstalt  fur  Luft-  und  Raumfahrt  e.V.  decided  in  1987  to 
develop  a  solar  thermal  desalination  project  at  the  Plataforma  Solar  de  Almeria.  In  July, 
1988,  the  first  start-up  of  a  solar  multi-effect  distillation  system  took  place  at  the 
Plataforma  Solar  de  Aimeria,  a  solar  research  centre  located  in  southeastern  Spain,  near 
Aimeria.  Its  evaluation  was  finished  in  December,  1990.  The  plant,  known  as  “So  1-14” 
was  built  and  connected  to  the  previous  existing  solar  facilities  as  a  result  of  a 
Spanish-German  project  [64].  The  performance  of  this  solar  desalination  system  was 
evaluated  by  Garcia-Rodriguez  and  Gomez-Camacho  [64]  using  energy  and  exergy 
analysis  method  in  order  to  propose  possible  improvements  of  the  system,  as  presented 
partly  below. 

In  the  exergetic  analysis  of  the  system  indicated  in  Fig.  15  [64],  Eq.  (15)  was  used  and  the 
results  obtained  on  a  component  basis  are  listed  in  Table  3  [64].  It  is  obvious  that  solar 
multi-effect  distillation  with  14  cells  had  the  lowest  exergetic  efficiency,  while  the  thermal 
storage  had  the  highest  one. 
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Fig.  15.  Solar  multi-effect  distillation  (MED)  system  at  Phase  I  of  the  solar  thermal  desalination  project  [64].  (1) 
Solar  input,  (2)  optic  and  thermal  losses,  (3-6)  synthetic  oils  at  various  temperatures,  (7)  thermal  losses,  (8) 
saturated  steam,  (9)  saturated  water,  (10)  fresh  water  (product),  (11)  thermal  losses,  (12)  seawater,  (13)  cooling 
outlet  stream,  (14)  blowdown  and  (15)  auxiliary  power  consumption. 


Table  3 

Exergetic  performance  of  the  equipment  at  Phase  I  of  the  solar  thermal  desalination  project  [64] 


Equipment 

F  (kW) 

P  (kW) 

F  —  P  (kW) 

£  =  P/F  (%) 

Solar  collector 

336.9 

58.63 

278.3 

17.4 

Thermal  storage 

58.63 

=  58.63 

=  0 

^  100 

Boiler 

58.63 

26.21 

32.42 

44.7 

Solar  multi-effect  distillation  (MED) 

34.91 

4.992 

29.92 

14.3 

Global  system 

345.6 

4.992 

340.6 

1.4 

4. 1.1. 7.  Solar  thermal  power  plants.  The  solar  thermal  power  system  in  general  can  be 
considered  as  consisting  of  two  subsystems,  namely,  the  collector-receiver  subsystem  and 
heat  engine  subsystem  as  illustrated  in  Fig.  16  [12].  The  collector-receiver  circuit  consists 
of  a  number  of  parabolic  trough  collectors  arranged  in  modules  operating  in  tracking 
mode.  The  main  components  of  the  heat  engine  circuit  are  a  heater  (here  the  boiler  heat 
exchanger),  a  turbine  having  two  stages,  a  condenser,  pump  and  a  regenerator.  Exergy 
analysis  is  conducted  in  two  main  subsystems  along  with  their  components,  as  given  by 
Singh  et  al.  [12]. 

The  overall  exergy  efficiency  of  the  collector-receiver  circuit  is  calculated  in  a  similar 
manner  given  by  Eq.  (22)  as  follows  [12]: 

_  Exu,c0l-rec  (r 

fco,col-rec  —  y-,  .  vjy  1 

Expar?s 

with  the  useful  exergy  of  the  whole  parabolic  collector 


Exu>col-rec  —  77par,colExu 

and  the  exergy  received  by  the  parabolic  collector  given  as: 


Expar,s  =  4(l-^\  (71) 

where  /7pa r,Coi  is  the  number  of  the  parabolic  collectors,  Qs  is  the  transferred  solar  heat  and 
Ts  is  the  solar  temperature.  The  exergy  received  by  the  parabolic  collector  may  also  be 
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Fig.  16.  Schematic  of  a  solar  power  thermal  plant  drawn  by  Singh  et  al.  [12].  (I)  Parabolic  through  collector,  (II) 
boiler  heat  exchanger,  (III)  high-pressure  turbine,  (IV)  low-pressure  turbine,  (V)  condenser,  (VI)  valve,  (VII) 
regenerator. 


calculated  using  Eqs.  (25)-(26),  and  the  relations  given  in  Table  1.  The  useful  exergy  of  the 
one  parabolic  collector  (Exu)  may  be  found  using  Eq.  (23b)  or  (23c). 

The  total  efficiency  of  the  heat  engine  subsystem  consisting  of  the  boiler  heat  exchanger 
and  the  heat  engine  cycle  is  calculated  as  follows: 


^net 

£h,eng  — 

-*^Xh,eng 

with  the  available  exergy  of  the  working  fluid  of  the  heat  engine  cycle  given  as: 


Exh>eng  =  QU(  1  -  -3®-) ,  (73) 

\  1  Ran,0  J 

where  Qu  is  the  useful  transferred  heat,  Wnet  is  the  net  work  done  by  the  heat  engine  cycle 
and  TRan  0  is  the  ambient  temperature  of  the  Rankine  cycle  [12]. 

The  overall  exergy  efficiency  of  the  solar  thermal  power  system  is  the  product  of  the  two 
separate  efficiencies  of  the  subsystems,  as  given  below  [12]: 


Ex 


par,s 


£h,eng  — 


(74) 
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You  and  Hu  [65]  investigated  the  optimal  thermal  and  exergetic  efficiencies  of  a 
medium-temperature  solar  thermal  power  system  using  energy  and  exergy  analysis 
method.  A  simple  reheat-regenerative  arrangement  is  schematically  shown  in  Fig.  17  [65]. 

The  exergetic  efficiency  of  the  cycle  (including  the  exergy  loss  in  the  heating  process)  is  as 
follows  [65]: 


_  fE turbine 

£cyde  - 


Exq  is  calculated  using  Eq.  (23c)  where  Tw  out  and  Twout  are  the  temperatures  of  the 
solar  heat  carrier  entering  and  exiting  the  heat  exchanger(s),  respectively. 

The  exergetic  efficiency  of  the  collector  is  calculated  from  [65]: 


£par,s  — 


Exq 

Exs 


with 


'  4  T 


am 


3  77 
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0.28/) 


The  overall  exergy  efficiency  of  the  combined  system  of  the  collector  and  the  power 
cycle  is  the  product  of  the  two  separate  efficiencies  of  the  subsystems,  as  given  by  the 
equation  [65]. 


£sys  —  £cycle£par,s-  tV*/ 

Singh  et  al.  [12]  calculated  the  irreversibilities  of  the  solar  thermal  power  plant  illustrated 
in  Fig.  16  [12].  It  was  found  that  the  collector  and  receiver  had  the  highest  irreversibilities 
and  the  efforts  to  reduce  exergy  loss  should  be  made  in  this  assembly.  The  component 
exergy  efficiency  values  were  as  follows:  Collector:  29.03%,  receiver:  68.38%,  collector- 
receiver:  19.85%,  boiler  heat  exchanger:  96.42%,  heat  engine:  66.57%  and  overall  system: 
12.74%. 

You  and  Hu  [65]  reported  based  their  analysis  that  the  reheat-regenerative  arrangement 
was  suitable  for  the  medium-temperature  solar  thermal  power  generation  purpose.  There 


5  (375°C)  i  (147.859kPa,  375°C) 


Fig.  17.  Schematic  of  one-stage  reheat-regenerative  power  cycle  and  the  parabolic  trough  collector  drawn  by  You 
and  Hu  [65]. 
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was  an  optimum  saturation  temperature  in  the  boiler  or  an  optimum  temperature  the  fluid 
entering  the  solar  field  (with  the  same  exiting  temperature)  at  which  both  the  thermal 
efficiency  and  the  exergetic  efficiency  of  the  combined  system  of  the  reheat-regenerative 
cycle  and  the  solar  collector  hit  the  maximum  values  simultaneously. 

4.1.2.  Photovoltaics 

Exergetic  evaluation  of  PVs  has  been  performed  by  some  investigators  [66,67]  in  hybrid 
(PV/thermal)  systems  as  a  part  of  the  system.  This  evaluation  is  explained  in  Section  4.1.3 
in  more  detail. 

Electrical  energy  is  not  affected  by  ambient  conditions  and  therefore  is  equivalent  in 
work.  If  global  irradiance  is  /,  energetic  efficiency  of  the  solar  cell  is  7/scell,  the  instantaneous 
electrical  exergy  is  then  as  follows  [66]. 

Exe  =  f/scell7  =  fiscell  I,  (79) 

where  esceii  is  the  exergetic  efficiency  of  the  solar  cell. 


4.1.3.  Hybrid  (PV/thermal)  solar  collectors 

The  exergy  efficiency  of  a  hybrid  solar  collector,  that  generates  both  electric  power  and 
heat,  may  be  calculated  as  follows  [67]. 


^PV/Thermal  — 


o 


conver 


/  +  Ex 


solar 


with 


Ex 


solar  — 


Tfluid  ~  Tp 

Tfluid 


Ssolar 


and 


Exi  =  0.95/,  (82) 

where  gconYQV  is  the  conversion  efficiency,  /  the  global  irradiance  (W/m2),  Exsoiar  the  exergy 

O  •  O  * 

of  heat  (W/m  ),  Exj  the  exergy  of  global  irradiance  (W/m  ),  gsolar  the  collected  solar  heat 
amount  per  unit  time  per  panel  area  (W /m  )  and  /fluid  the  supply  temperature  of  collector 
fluid  (K). 

Fujisawa  and  Tani  [66]  defined  the  synthetic  exergy  of  the  PV/T  collector  as  the  total 
value  of  the  electrical  and  thermal  exergies  as  follows: 

EXpy/T  =  Exe  T  EXiflermal  —  (^e  T  ^thermal)  ~  ^PV/T I  ($3) 

with 


Exe  =  sj  (84) 

and 

^  —  ^thermalE  ($5) 

Fujisawa  and  Tani  [66]  designed  and  constructed  a  PV-thermal  hybrid  collector  on  their 
campus.  The  collector  consisted  of  a  liquid  heating  fiat-plate  solar  collector  with  mono-Si  PV 
cells  on  substrate  of  non-selective  aluminum  absorber  plate,  The  collector  area  was 
1.3  x  0.5  m2.  From  the  annual  experimental  evaluation  based  on  exergy,  Fujisawa  and  Tani 
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concluded  that  the  PV/T  collector  could  produce  higher  output  density  than  a  unit  PV  module 
or  liquid  heating  fiat-plate  solar  collector.  Using  exergetic  evaluation,  the  best  performance  of 
available  energy  was  found  to  be  that  of  the  coverless  PV/T  collector  at  80.8kWh/yr,  the 
second  to  be  the  PV  module  at  72.6kWh/yr,  the  third  to  be  the  single-covered  PV/T  collector 
at  71.5kWh/yr,  and  the  worst  to  be  the  flat-plate  collector  at  6.0kWh/yr. 

Saitoh  et  al.  [67]  performed  field  measurements  from  November  1998  to  October  1999  at 
a  low  energy  house  at  Hokkaido  University.  A  system  diagram  of  the  experimental 
equipment  is  illustrated  in  Fig.  18  [67].  The  power  generated  was  measured  by  giving  the 
variable  load  for  the  maximum  power  point  tracking.  Two  hybrid  solar  panels  were 
connected  in  parallel,  while  the  brine  was  supplied  at  constant  temperature  by  a  fluid 
supply  system  with  a  circulating  pump.  The  volumetric  flow  rate  per  panel  was  fixed  at  1 1/ 
min.  The  tilt  angle  of  the  panel  was  30°,  which  gave  the  maximum  annual  global 
irradiance.  Except  for  winter,  the  mean  conversion  efficiency  of  array  and  the  collector 
efficiency  were  stable  at  8-9%  and  25-28%,  respectively.  The  dependency  on  solar  energy 
was  46.3%.  Energy  and  exergy  efficiency  values  of  single-junction  crystalline  silicon  PV  are 
illustrated  in  Table  4  [67].  The  exergy  efficiencies  of  solar  collector,  PV  and  hybrid  solar 
collector  were  found  to  be  4.4%,  11.2%  and  13.3%,  respectively.  Solar  collector  had  the 
lowest  value.  In  terms  of  energy  quality,  the  hybrid  solar  collector  was  very  effective  [67]. 

4.2.  Wind  energy  systems 

The  exergy  of  fluid  currents  is  entirely  due  to  the  kinetic  energy  of  the  fluid.  The 
maximum  work  is  extracted  from  a  moving  fluid  when  the  velocity  is  brought  to  zero 
relative  to  the  reference  state.  The  kinetic  energy  of  the  fluid  and  its  exergy  have  the  same 
numerical  value  [14]. 

'/'fluid  =  \pV0’  (86) 

where  p  is  the  density  of  the  wind  and  V0  is  the  wind  speed. 


Fig.  18.  System  diagram  of  experimental  equipment  drawn  by  Saitoh  et  al.  [67]. 
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Table  4 

Comparison  of  energy  and  exergy  efficiency  values  of  solar  collector,  photovoltaic  and  hybrid  collector  [67] 


Type  of  efficiency 

Solar  collector 

Photovoltaic 

Hybrid  collector 

Energy  efficiency  (g,  %) 

Heat 

46.2 

— 

32.0 

Power 

— 

10.7 

10.6 

Total 

46.2 

10.7 

42.6 

Exergy  efficiency  (e,  %) 

Heat 

4.4 

— 

2.1 

Power 

— 

11.2 

11.2 

Total 

4.4 

11.2 

13.3 

Koroneos  et  al.  [11]  presented  two  diagrams  where  utilization  of  the  wind’s  potential  in 
relation  to  the  wind  speed  and  exergy  losses  in  the  different  components  of  a  wind  turbine 
(i.e.,  rotor,  gearbox  and  generator)  were  included.  They  concluded  as  follows:  “According 
to  Betz’s  law,  wind  turbine  can  take  advantage  of  up  to  60%  of  the  power  of  the  wind. 
Nevertheless,  in  practice,  their  efficiency  is  about  40%  for  quite  high  wind  speeds.  The  rest 
of  the  energy  density  of  the  wind  not  obtainable  is  exergy  loss.  This  exergy  loss  appears 
mainly  as  heat.  It  is  attributed  to  the  friction  between  the  rotor  shaft  and  the  bearings,  the 
heat  that  the  cooling  fluid  abducts  from  the  gearbox,  the  heat  that  the  cooling  fluid  of  the 
generator  abducts  from  it  and  the  thyristors,  which  assist  in  smooth  starting  of  the  turbine 
and  which  lose  1-2%  of  the  energy  that  passes  through  them.” 

As  for  as  studies  conducted  on  exergetic  assessment  of  wind  energy  systems,  Sahin  et  al. 
[68]  developed  a  new  exergy  formulation  for  wind  energy,  which  was  more  realistic  and 
also  accounted  for  the  thermodynamic  quantities  enthalpy  and  entropy.  The  relation 
developed  for  the  total  exergy  for  wind  energy  is  as  follows,  while  the  input  and  output 
variables  for  the  system  considered  are  illustrated  in  Fig.  19  [68]. 

Ex  =  £gen  +  mCp(T2  -  Tx)  +  mTat  (cp\n^-R  In  §  -  %*)  (87) 

\  3  I  E\  1  at  J 


with 


aoss  =  mCp(rat  -  rave),  (88) 

where  Egen  is  the  generated  electricity  (i.e.,  the  total  kinetic  energy  difference),  Tat  is  the 
atmospheric  temperature,  P2  is  the  pressure  at  the  exit  of  the  wind  turbine  for  a  wind  speed 
V2  and  Pi  is  the  pressure  at  the  inlet  of  the  wind  turbine  for  a  wind  speed  F1?  gloss 
represents  heat  losses  of  wind  turbine  and  Tave  is  the  mean  value  of  input  and  output  wind 
chill  temperatures. 

Sahin  et  al.  [68]  reported  as  follows:  “The  exergy  analysis  of  wind  energy  shows  that 
there  are  significant  differences  between  energy  and  exergy  analysis  results.  According  to 
one  classical  wind  energy  efficiency  analysis  technique,  which  examines  capacity  factor,  the 
resultant  wind  energy  efficiency  is  overestimated  relative  to  what  it  really  is.  The  capacity 
factor  normally  refers  to  the  percentage  of  nominal  power  that  the  wind  turbine  generates. 
Average  differences  between  energy  and  exergy  efficiencies  are  approximately  40%  at  low 
wind  speeds  and  up  to  approximately  55%  at  high  wind  speeds”. 
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Fig.  19.  Wind  turbine  and  representative  wind  energy  input  and  output  variables  drawn  by  Sahin  et  al.  [68]. 


Ozgener  and  Ozgener  [69]  performed  exergy  analysis  of  a  wind  turbine  system  (1.5  kW) 
located  in  Solar  Energy  Institute  of  Ege  University  (latitude  38.24N,  longitude  27.50E), 
Izmir,  Turkey.  They  reported  that  exergy  efficiency  changed  between  0%  and  48.7%  at 
different  wind  speeds  based  on  a  dead  state  temperature  of  25  °C  and  a  atmospheric 
pressure  of  101.325  kPa,  considering  pressure  differences  between  state  points.  Depending 
on  temperature  differences  between  state  points,  exergy  efficiencies  were  found  to  be 
0-89%. 

The  exergetic  efficiency  of  a  turbine  is  defined  as  a  measure  of  how  well  the  stream 
exergy  of  the  fluid  is  converted  into  useful  turbine  work  output  or  inverter  work  output.  By 
using  inverter  work  output  as  useful  work  and  applying  this  to  wind  turbine,  the  following 
relations  are  obtained  [69]. 


8 


W 


e,io 


Exi  -  Ex2 


with 


Exdest  =  (Exi  -  Ex2)  -We,  (90) 

where  WQ-l0  is  the  power  at  inverter  output  and  Eldest  is  the  exergy  destruction,  while  Exi 
and  Ex2  are  exergy  rates  calculated  using  Eqs.  (12)  and  (13). 


4.3.  Geothermal  energy  systems 

The  specific  exergy  of  a  geothermal  site  depends  on  the  temperature  at  the  collector,  the 
pressure  and  composition  of  fluid  in  the  reservoir,  and  the  ability  of  the  surrounding  rock 
or  fluid  to  transport  thermal  energy  to  the  collector.  Utilization  of  a  geothermal  site  often 
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changes  the  thermal  or  fluid  properties  of  the  reservoir  to  a  new  equilibrium  dependent  on 
the  method  of  energy  extraction  [14]. 


4.3.1.  Classification  of  geothermal  resources  by  exergy 

Geothermal  resources  have  been  classified  as  low,  intermediate  or  high  enthalpy 
resources  according  to  their  reservoir  temperatures.  The  temperature  ranges  used  for  these 
classifications  are  arbitrary  and  they  are  not  generally  agreed  upon.  Temperature  is  used  as 
the  classification  parameter  because  it  is  the  earliest  to  measure  and  understand.  In 
addition,  temperature  or  enthalpy  alone  can  be  ambiguous  defining  a  geothermal  resource 
because  two  independent  thermodynamic  properties  are  required  to  define  the  thermo¬ 
dynamic  state  of  fluid.  Geothermal  energy  is  already  in  the  form  of  heat,  and  from  the 
thermodynamic  point  of  view,  work  is  more  useful  than  heat  because  not  all  the  heat  can 
be  converted  to  work.  Therefore,  geothermal  resources  should  be  classified  by  their  exergy, 
a  measure  of  their  ability  to  do  work  [70,71]. 

Lee  [71]  proposed  a  new  parameter,  the  so-called  specific  exergy  index  (SExI)  for  better 
classification  and  evaluation  as  follows: 


SExI  = 


^brine  273.1 6$brine 

U92 


which  is  a  straight  line  on  an  h—s  plot  of  the  Mollier  diagram.  Straight  lines  of  SExI  =  0.5 
and  SExI  =  0.05  can  therefore  be  drawn  in  this  diagram  and  used  as  a  map  for  classifying 
geothermal  resources  by  taking  into  account  the  following  criteria: 


•  SExI  <0.05  for  low-quality  geothermal  resources, 

•  0.05^  SExI  <0.5  for  medium-quality  geothermal  resources,  and 

•  SExI  ^0.5  for  high-quality  geothermal  resources. 


Here,  the  demarcation  limits  for  these  indices  are  exergies  of  saturated  water  and  dry 
saturated  steam  at  1  bar  absolute. 

In  order  to  map  any  geothermal  field  on  the  Mollier  diagram  as  well  as  to  determine  the 
energy  and  exergy  values  of  the  geothermal  brine,  the  average  values  for  the  enthalpy  and 
entropy  are  then  calculated  from  the  following  equations  [72]: 


^brine  — 


V„  riiw,ihwj 

En  • 

/  mwJ 


^brine  — 


En  • 

/  mw,iS 


w,i 


En  • 

m^i 


(92) 

(93) 


The  above  classification  has  been  used  by  some  investigators  in  their  studies.  In  this 
regard,  Quijano  [72]  performed  exergy  analysis  of  the  Ahuachapan  and  Berlin  geothermal 
fields.  Plotting  the  thermodynamic  conditions  of  both  fields  on  the  Molliere  diagram  and 
computing  the  specific  exergy  index  (SEI),  the  fields  were  classified  on  the  medium  high 
exergy  zone.  Ozgener  et  al.  [73]  and  Baba  et  al.  [74]  applied  the  SExI  to  the  Balcova 
geothermal  field,  which  is  located  in  the  western  part  of  Turkey,  and  is  endowed  with 
considerably  rich  geothermal  resources.  The  number  of  the  wells  in  operation  in  this  field 
may  vary  depending  on  the  heating  days  and  operating  strategy.  Taking  into  account  the 
eight  productive  wells,  which  were  in  operation  at  the  date  considered  and  using  Eq.  (91) 
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together  with  the  values  given  in  Table  5  [74],  the  SExI  was  found  to  be  0.07.  This 
represented  that  the  field  fell  into  the  medium  quality  geothermal  resources  according  to 
the  classification  of  Lee  [71].  For  classifying  the  Balcova  geothermal  field  according  to  its 
reservoir,  the  weighted  average  value  of  the  temperature  (rbrine)  was  obtained  to  be 
1 13.9  °C  by  using  the  similar  relation  given  in  Eq.  (91).  This  field  fell  into  the  intermediate 
enthalpy  resources  according  to  the  classifications  of  Muffler  and  Cataldi  [75]  and 
Benderitter  and  Cormy  [76],  while  it  is  characterized  as  low  enthalpy  resources  according 
to  the  classification  of  Hochstein  [77].  This  clearly  indicated  that  the  classification  with 
reference  to  SExI  was  more  meaningful  as  there  was  no  general  agreement  on  the  arbitrary 
temperature  ranges  used  in  the  classification  of  geothermal  resources  by  temperature. 
Ozgener  et  al.  [78,79]  also  determined  the  SExI  values  for  the  Gonen  and  Salihli 
geothermal  fields  in  Turkey  as  0.025  and  0.049,  respectively.  They  concluded  that  both 
these  fields  fell  into  the  low  quality  geothermal  resources  according  to  the  classification  of 
Lee.  Etemoglu  and  Can  [80]  classified  geothermal  resources  in  Turkey  based  on  SExI.  In 
their  analysis,  the  triple  point  of  water  was  selected  as  the  dead-state  conditions  since  the 
enthalpy  and  entropy  values  were  zero.  SExI  values  of  geothermal  resources  over  0.56  were 
classified  as  high  exergy  resources,  SER<  0.05  low  exergy  and  between  of  these  values  were 
medium  exergy  resources. 

4.3.2.  Direct  utilization 

4.3.2. 1.  District  heating.  Although  various  studies  have  been  undertaken  to  perform 
energy  and  exergy  analyses  of  geothermal  systems,  the  exergy  analysis  of  geothermal 
district  heating  systems  (GDHSs)  based  on  the  actual  data  was  performed  by  Ozgener 
et  al.  [81]  for  the  first  time  (to  the  best  of  authors’  knowledge).  In  this  context,  the  first 
study  included  a  comprehensive  energy  and  exergy  analysis  conducted  on  the  Balcova 
geothermal  district  heating  system  located  in  the  town  Balcova,  Izmir,  Turkey,  using  actual 
data  sets  obtained  as  of  January  1,  2003  for  system  analysis,  performance  evaluation  and 
efficiency  improvement.  Since  then,  Ozgener  et  al.  have  conducted  various  studies  on 
exergetic  [82,83]  and  exergoeconomic  analyses  [84,85]  of  GDHSs,  namely  Balcova,  Salihli 
and  Gonen  GDHSs  installed  in  the  provinces  of  Izmir,  Manisa  and  Balikesir,  Turkey 


Table  5 

Some  thermodynamic  properties  of  the  wells  in  the  Balcova  geothermal  field  in  Izmir,  Turkey  (as  of  1  January 
2003)  [74] 


Code  of  well 

Temperature  T 
(°C) 

Mass  flow  rate  m 
(kg/s) 

Pressure  P 
(kPa) 

Specific  enthalpy 
h  (kJ/kg) 

Specific  entropy 

5  (kJ/kg  K) 

BD2 

131.6 

18.31 

413.4 

552.98 

1.6507 

BD3 

120 

18.65 

522.83 

503.72 

1.5270 

BD4 

135 

25.32 

539.04 

567.57 

1.6863 

BD5 

111.4 

3.74 

607.95 

467.35 

1.4333 

BD7 

110.6 

20.01 

496.49 

463.89 

1.4245 

B4 

98.8 

13.16 

368.81 

413.99 

1.2926 

B5 

103.6 

23.10 

332.34 

434.20 

1.3467 

B10 

96 

25.90 

364.76 

402.19 

1.2607 

BD8 

(reinjection) 

66.1 

112.02 

474.20 

276.8 

0.9061 

Average  values 

113.9 

148.19 

— 

477.87 

1.4578 
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as  well  as  on  investigating  the  effect  of  dead  state  temperature  on  the  energy  and 
exergy  efficiencies  of  GDHSs  [86,87].  They  also  performed  a  key  review  on  the 
performance  improvement  aspects  of  GDHSs  [88]  and  prepared  a  general  paper  on 
geothermal  energy  including  the  exergetic  evaluation  of  GDHSs  for  the  Encyclopedia  of 
Energy  Engineering  [89]. 

The  following  relations  have  been  summarized  from  the  previous  studies  performed  by 
Ozgener  et  al.  [81-91]  as  follows: 

A  schematic  diagram  of  the  simplified  GDHS,  where  user  system  is  heated  by 
geothermal  energy,  is  illustrated  in  Fig.  20  [88].  As  can  be  seen  in  this  figure,  a  GDHS 
consists  mainly  of  three  circuits,  namely  (a)  energy  production  circuit  (EPC;  geothermal 
well  loop),  (b)  energy  distribution  circuit  (EDC;  district  heating  distribution  network,  and 
(c)  energy  consumption  circuit  (ECC).  In  evaluating  the  system  performance  from  the 
exergetic  point  of  view,  the  two  items  of  the  whole  system  play  a  big  role.  These  include:  (i) 
the  geothermal  brine  exergy  inputs  from  the  production  field  and  (ii)  exergy  destructions  in 
the  system  taking  place  as  the  exergy  of  the  fluid  lost  in  the  heat  exchanger,  the  natural 
direct  discharge  of  the  system  (pipeline  losses),  and  the  pumps. 

For  a  general  steady-state,  steady-flow  process,  the  two  balance  equations,  namely  mass 
and  exergy  balance  equations,  may  be  written  as  follows: 

Using  Eq.  (1),  the  mass  balance  equations  for  two  various  GDHSs,  namely  GDHSs  with 
reinjected  thermal  water  (i.e.,  Balcova  GDHS)  and  without  reinjected  thermal  water 
(i.e.,  Salihli  GDHS),  may  be  written  as  follows,  respectively  [88]: 

n 

T  "VT  -  mr  -  m&  =  0  (94) 

i=\ 

and 


n 

^w,T  -  rhd  =  o,  (95) 

i—  1 


EPC 


(a)  Energy  production 
cycle  (EPC) 


EPC 


U 

Ph 

w 


Central  Heating 
Station 


EDC 


(b)  Energy  distribution 
cycle  (EDC) 


EDC _ ^  EDC 

Circulating  pump 


ECC 

(c)  Building  substations - 

(Energy  consumption 

circuit  -ECC-)  - 

^ECC  ^ _ ECC 

Circulating  pump 


Reinjection  well 

Fig.  20.  Basic  simplified  geothermal  district  heating  energy  system  schema  [90]. 
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where  mwj  is  the  total  mass  flow  rate  at  wellhead,  mr  is  the  flow  rate  of  the  reinjected 
thermal  water  and  m d  is  the  mass  flow  rate  of  the  natural  direct  discharge. 

The  exergy  destructions  in  the  heat  exchanger,  pump  and  the  system  itself  are  calculated 
as  follows  [88]: 

Ex^es^j-iE  =  EXin  Exout  =  EXdest,  (96) 


and 


EXdest,pump  —  fE  pump  (Exout  EXin) 


EXdest  , system  —  E  EXdest,HE  T  ^  ^  EXdest, 


(97) 


pump  • 


(98) 


Using  Eq.  (15),  the  exergy  efficiencies  of  two  various  GDHSs  with  reinjected  thermal 
water  and  without  reinjected  thermal  water  are  found  from  the  following  relations, 
respectively  [88]. 


£GDHS,wr  = 


Exusefui5HE 

EXbrine 


=  1  - 


EXdest, system  T  Exrejnjected  +  Exnaturai  discharged 

EXbrine 


(99) 


Exusefui]dE  1  EXdest  ,  system  +  Exnaturai  discharged 

%DHS,wor  =  - • -  =  1 - • - 


EXbrine 


Ex^rine 


(100) 


where  Eldest, system  is  the  exergy  destructions  in  the  system  taking  place  as  the  exergy  of  the 
fluid  lost  in  the  heat  exchanger  and  the  pumps,  Exreinjected  is  the  exergy  loss  due  to  the 
reinjected  water  and  Exnaturai  discharged  is  the  exergy  destruction  due  to  the  natural  direct 
discharge  of  the  system  (due  to  a  significant  amount  of  water  leaks;  pipeline  losses). 

Table  6  [81-83,86,87,90,91]  shows  a  comparison  of  exergy  efficiency  values  for  GDHSs 
installed  in  Turkey.  As  can  be  seen  in  this  table,  exergy  efficiencies  range  from  42.94%  to 
66.13%  depending  on  the  various  dead  state  temperatures.  As  expected,  the  lower  the 
reference  state  temperature,  the  significantly  larger  the  exergy  losses  in  the  system  pipeline 
and  heat  exchangers.  However,  exergy  rates  of  pipeline  losses  of  the  system,  heat 
exchangers,  and  useful  exergy  increase  considerably.  The  reason  for  this  rapid  rise  in 
exergy  rates  is  due  to  a  decrease  in  the  ambient  temperature  [88]. 

Ozgener  et  al.  [88]  also  studied  on  the  effect  of  varying  dead  state  temperatures  on  the 
exergy  efficiencies  of  the  whole  GDHS.  Exergy  efficiency  values  were  correlated  with  high 
correlation  coefficients  for  the  Balcova,  Gonen  and  Salihli  GDHSs.  Of  these,  the 
correlation  for  the  Balcova  GDHS  is  given  below  [87] 


^Balcova,  gdhs  =  0.000006  —  0.000057^m  —  0.1279Tam  +  45.21,  (101) 


where  Tam  is  the  ambient  temperature  in  °C. 


4. 3. 2. 2.  Ground-source  (geothermal)  heat  pumps.  Ground-source  heat  pumps  (GSHPs) 
have  been  used  for  years  in  the  developed  countries  due  to  their  higher  energy  utilization 
efficiencies  compared  with  conventional  heating  and  cooling  systems.  There  are  two  main 
types  of  GSHP  system,  (a)  ground-coupled  (vertical  or  horizontal)-closed  loop  and  (b) 
water  source-open  loop.  Ground-coupled  heat  pumps  are  known  by  a  variety  of  names. 
These  include  GSHPs,  earth-coupled  heat  pumps,  earth  energy  heat  pumping  systems, 
earth  energy  systems,  ground-source  systems,  geothermal  heat  pumps,  closed-loop  heat 


Table  6 

Comparison  of  various  geothermal  district  heating  systems  in  terms  of  exergy  destructions  and  exergy  efficiencies  [8 1—83,90,9 1 ,86,87] 


Name  of 

GDHS 

Location/ 

country 

Dead  state 

temperature 

(°C) 

Date  of  data 
used 

Total  exergy 
input  (kW) 

Exergy  destructions/losses  (°/c 
input) 

i  of  the  total  (brine)  exergy 

Exergy 

efficiency 

■  (%) 

References 

Pumps 

Heat 

exchangers 

Thermal 

reinjection 

Natural 
direct 
discharge 
(pipe  line) 

Balcova 

GDHS 

Izmir /Turkey 

13.1 

1  January 
2003 

9164.29 

3.22 

7.24 

22.66 

24.1 

42.94 

[81] 

Gonen 

GDHS 

Balikesir/ 

Turkey 

6 

1  February 
2004 

2657.5 

14.81 

7.11 

12.96 

1.06 

64.06 

[82] 

Salihli 

GDHS 

Manisa/ 

Turkey 

2.9 

1  February 
2004 

2564 

2.22 

17.90 

— 

20.44 

59.44 

[83] 

Balcova 

GDHS 

Izmir /Turkey 

13. la 

2  January 
2004 

14808.15 

1.64 

8.57 

14.84 

28.96 

46.00 

[84] 

Gonen 

GDHS 

Balikesir/ 

Turkey 

10b 

1  February 
2004 

2333.33 

17.45 

8.19 

7.37 

0.86 

66.13 

[85] 

Salihli 

GDHS 

Manisa/ 

Turkey 

25c 

1  February 
2004 

1524 

4.27 

31.17 

— 

8.98 

55.58 

[88] 

Balcova 

GDHS 

Izmir/Turkey 

11. 4d 

2  January 
2004 

14390 

1.74 

8.83 

14.18 

28.7 

46.55 

[89] 

aThe  study  was  conducted  for  four  various  dates  ,  namely  1  January  2003 , 2  December  2003,  2  January  2004  and  2  February  2004  at  dead  state  temperatures  of 
13.1,  12,  10.4  and  6.5  °C,  respectively. 

bA  parametric  study  was  performed  o  investigate  how  varying  reference  how  varying  reference  temperature  from  0  to  20  °C  would  affect  the  exergy  efficiencies. 
cThe  study  was  conducted  for  1  December  2003,  1  January  2004  and  1  February  2004  at  dead  state  temperatures  of  10.8,  9.7  and  2.9,  respectively.  In  addition,  two 
parametric  expressions  of  exergy  efficiencies  were  developed  as  a  function  of  the  reference  temperature,  by  varying  the  dead  state  temperature  from  0  to  25  °C. 

dThe  reference-state  temperature  value  of  1 1.4  °C  was  obtained  from  the  average  ambient  temperature  data  measured  between  2000  and  2004  for  second  January  in 
Izmir.  The  dead  state  temperature  was  considered  varying  from  0  to  25  °C  to  conduct  a  parametric  study. 
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pumps,  solar  energy  heat  pumps,  geoexchange  systems,  geosource  heat  pumps,  and  a  few 
other  variations  [92]. 

Fig.  21  [93]  illustrates  a  schematic  diagram  of  a  GSHP  system  in  the  heating  mode  as  an 
illustrative  example.  The  system  was  designed  and  installed  in  the  Solar  Energy  Institute  of 
Ege  University,  Izmir,  Turkey.  The  system  was  commissioned  in  May  2000,  while  GSHP 
systems  have  been  applied  to  the  Turkish  residential  buildings  since  1997. 

The  performance  tests  have  also  been  performed  since  then.  This  system  has  been 
recently  integrated  with  a  solar  collector  and  a  greenhouse  [94].  It  mainly  consists  of  three 
separate  circuits:  (i)  the  ground  coupling  circuit  (brine  circuit  or  water-antifreeze  solution 
circuit),  (ii)  the  refrigerant  circuit  (or  a  reversible  vapor  compression  cycle)  and  (iii)  the  fan 
coil  circuit  (water  circuit). 

In  Table  7  [95],  Eqs.  (1 02)— (1 26)  used  in  modeling  and  analyzing  GSHP  systems  are 
summarized  [95].  These  relations  have  been  obtained  for  each  of  the  GSHP  components 
illustrated  in  Fig.  21  using  mass,  energy,  entropy  and  exergy  balance  equations  as  well  as 
the  exergy  destructions  obtained  using  exergy  balance  equations  given  in  Section  3). 


Ground  level 


Fig.  21.  The  main  components  and  schematic  of  a  ground-source  heat  pump  system  [93].  (I)  Compressor, 
(II)  condenser,  (III)  capillary  tube,  (IV)  evaporator,  (V)  fan-coil  unit,  (VI)  ground  heat  exchanger. 


Table  7 

Mass,  energy,  exergy  and  entropy  balance  relations  of  a  GSHP  system  [95] 


No  Element 


Mass  analysis  using  Eq.  (1)  (the  The  irreversibility  or  the  exergy  destroyed  using  Eq.  The  irreversibility  or  the  exergy  destroyed  using 

conversion  of  mass  principle)  (5a)  (the  exergy  balance)  entropy  balance  and  Eq.  (10)  (the  decrease  of  exergy 

principle) 


I 


II 


III 


IV 


V 


VI 


Compressor  777 1  —  ,n2  —  777  r 


Condenser  777 1  —  7772  —  777  r 

rhs  =  the  =  mw 

Throttling  valve  ^3  =  7774  =  777  r 

(capillary  tube) 


Evaporator 


m4  =  m\  =  mT 


Fan  coil  units 


777air,in  —  777  air,  out  —  777  air 


Ground  heat  777  7  —  777  8  —  777wa 

exchanger 


(102)  Exi  +  W  comp,e  =  Ex?a  +  7[ 

7 1  =  777r(lAl  </%)  T  fE comp 

(106)  Ex2a  +  Ex6  =  Ex3  +  Ex5  +  7n 

(107)  7n  =  mr(xl/2a  -  b3)  +  ^wOAe  “  •As) 

(111)  Ex3=Ex4+7ih 

7m  =  mr(iJ/3  -  i//4) 

(115)  Ex4  +  Ex8  =  Exi  +  Ex7  +  7iV 

7iv  =  /nr( <A4  -  <Al)  +  777wa(<As  - 


(103)  7t  =  r0mr(5i  -  52a) 

(104) 

(108)  7n  =  T0[mrCs'3  -  52a)  +  fhw(s5  -  j6)1 

(109) 

(112)  7m  =  T  othr(s4  —  s3) 

(113) 

(116)  7iv  =  T0[mr(s i  -  j4)  +  mwa0 7  -  s8)] 

(117) 


(119)  e  *  ,  fTmAil-T0\  ,  , 

Ex5  =  Ex6  +  Qfc  — ^ - )  +  /v 


T 


in, air 


(120)  7v  =  To 


777  w(v  -  55)  + 


2fc 


T 


in, air, 


7y  =  mw(lA5  -  <Ae)  -  0fc  1  - 


Tn 


(121) 


(123) 


TSoii  —  T  0 


\  =  Ex8  +  7  VI 
/ 

(124) 

< 

II 

O 

777wa(V  V)  T8 

7  soil 

7VI  =  Wwa(»A7  -  •As)  +  1 


Tn 


Tsoil 


(125) 


(105) 


(110) 


(114) 


(118) 


(122) 


(126) 
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For  further  details,  some  relevant  studies  published  should  be  taken  into  account 
[93-97]. 

Using  Eq.  (15),  exergy  efficiency  relations  for  the  GSHP  unit  and  the  whole  system  are 
written  as  follows: 


£o,r,hp  = 


EXheat  Exconcyn  ExconcjOU|; 


W 


a,m 


w 


(127) 


a,m 


£o,R,sys  — 


EXheat 


*Ea,in  +  W 


pumps 


+  rv. 


(128) 


an-coil 


COP 


EBP  = 


a,h 


COPc>h 


(129) 


Pt,hp 
£hp  =  — 

T,HP 


(130) 


^sys  — 


Pi. 


sys 


r  T 


(131) 


sys 


Exergy  analysis  has  been  applied  to  various  types  of  heat  pumps  by  a  number  of 
researchers  [i.e.,  41,50,51,56,98-100].  However,  the  studies  on  exergy  analysis  of  GSHPs  are 
relatively  few  in  numbers  [93-97,101-103].  Some  of  them  [i.e.,  43,104]  have  been  included  in 
the  exergoeconomic  (a  combination  of  exergy  and  economics)  analysis  of  GSHP  systems.  The 
first  detailed  exergetic  evaluation  of  a  GSHP  system  was  made  by  Hepbasli  and  Akdemir  [95] 
to  the  current  knowledge  of  the  author,  while  Ozgener  and  Hepbasli  [105]  also  conducted  a 
comprehensive  review  on  the  energy  and  exergy  analysis  of  solar  assisted  heat  pump  systems. 
In  this  regard,  the  performance  of  this  GSHP  system  with  a  50  m  vertical  32  mm  diameter  U- 
bend  ground  heat  exchanger  was  assessed  using  energy  and  exergy  analysis  methods.  Fig.  22 
[94]  illustrates  a  graphical  representation  of  the  exergy  balance,  namely  the  Grassmann 
diagram,  for  the  GSHP  system.  This  diagram,  on  the  other  hand,  gives  quantitative 
information  regarding  the  proportion  of  the  exergy  input  to  the  system  that  is  dissipated  in 
each  of  the  components.  As  can  be  seen  in  Fig.  22,  the  highest  irreversibility  occurred  in  sub- 
region  I,  the  motor-compressor  subassembly.  The  total  magnitude  of  the  losses  was  over  56% 
of  the  actual  power  input,  while  the  mechanical-electrical  losses  accounted  for  46%  of  that. 
These  losses  were  due  to  the  electrical,  mechanical  and  isentropic  efficiencies  and  emphasize 
the  need  for  paying  close  attention  to  the  selection  of  this  type  of  equipment,  since 
components  of  inferior  performance  could  considerably  reduce  the  overall  performance  of  the 
system.  The  second  largest  irreversibility  was  due  to  the  condenser.  This  was  partly  due  to  the 
large  degree  of  superheat  achieved  at  the  end  of  the  compression  process,  leading  to  large 
temperature  differences  associated  with  the  initial  phase  of  heat  transfer.  The  third  highest 
irreversibility  was  in  the  capillary  tube  due  to  the  pressure  drop  of  the  refrigerant  passing 
through  it.  Besides  this,  the  evaporator  had  the  lowest  irreversibility  on  the  basis  of  the  heat 
pump  cycle  [95]. 

Table  8  [93-97,103,106,107]  lists  various  ground  source  (geothermal)  heat  pumps 
analyzed  using  exergy  analysis  method  Five  types  of  exergy  efficiencies  given  in 
Eqs.  (127)— (1 3 1)  were  used  in  the  calculations.  It  is  obvious  from  this  table  that  the 
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Fig.  22.  Exergy  diagram  (the  Grassmann  diagram)  for  the  GSHP  system  [94].  (I)  Compressor,  (II)  condenser, 
(III)  capillary  tube,  (IV)  evaporator,  (V)  fan-coil  units,  (VI)  ground  heat  exchanger. 


exergy  efficiency  values  of  GSHP  units  vary  between  about  67%  and  92%  on  the  product/ 
fuel  basis  depending  on  the  various  dead  state  temperatures.  The  values  for  the  exergy 
improvement  rate  and  thermodynamic  parameters  calculated  using  Eq.  (16)  and 
Eqs.  ( 1 7)— (20)  were  also  included  in  this  table. 

4. 3. 2. 3.  Greenhouses.  Good  plant-growth  conditions  can  be  achieved  by  using  green¬ 
houses.  A  greenhouse  can  be  managed  to  protect  the  plants  by  creating  a  favorable 
environment,  allowing  intensive  use  of  soil,  and  helping  sanitary  plant  control.  A  variety  of 
heating  systems  are  being  used  in  present  day  greenhouses  to  meet  the  heating  and  cooling 
requirements.  As  an  alternative,  there  is  great  potential  in  employing  a  heat  pump  system 
for  greenhouse  air-conditioning  based  on  its  ability  to  perform  the  multi-function  of 
heating,  cooling  and  dehumidification  [11]. 

Although  various  studies  were  undertaken  to  evaluate  the  performance  of  solar  assisted 
heat  pump  systems,  to  the  best  of  author’s  knowledge,  two  studies  [1 1,14]  have  appeared  in 
the  open  literature  on  the  experimental  performance  testing  of  a  solar  assisted  vertical 
GSHP  system  for  greenhouse  heating.  The  first  one  [11]  was  related  to  an  exergetic 
assessment  at  a  fixed  dead  state  temperature,  while  the  second  one  [14]  included  a 
parametric  study  investigating  the  effect  of  varying  dead  state  temperature  on  the  exergy 
efficiencies. 


Table  8 

Various  ground  source  (geothermal)  heat  pumps  analyzed  using  exergy  analysis  method  [93-97,103,106,107] 


Location 

Type  of 

Capacity 

Dead 

Date 

Exergy  destructions  (rate  of  improvement  potential)  (%  of  the  whole  heat  pump  system) 

Eq.  no. 

Exergy 

Thermodynamic 

of  the 

heat  pump/ 

in  the 

state 

of  data 

for 

efficiency 

parameters  for 

system 

operating 

mode 

temp. 

used 

Ground  source  heat  pump  system  (Whole  system)1' 

exergy 

(%) 

compressor  (% 

built/ 

mode 

operated 

(°C) 

eff. 

of  the  whole  system) 

country 

(kW) 

Ground  source  heat  pump  unit  Other  components  Heat  Overall  %  5  £  / 


pump 

I  Com.  II  Con.  Ill  Thr.  IV  V  Fan-coil/  VI  GHE  VII  Cir.  VIII 

valve  Evap.  Con.  WT  /PHEx  pumps  Sol. 

col. 


Izmir  / 
Turkey 

Vertical 
GSHP/heating 
(for  space 
heating) 

3.41 

20 

Feb.  2001 

16.54 

27.91 

21.62 

7.83 

Izmir  / 
Turkey 

Vertical 
GSHP/heating 
(for  space 
heating) 

4.27 

2.2 

7  Jan. 2004 

51.83 

(57.12) 

8.00 

(3.03) 

8.36 

(2.43) 

9.42 

(14.30) 

Izmir/ 

Turkey 

Solar  assisted 
vertical  GSHP/ 
heating 

(for  greenhouse 
heating) 

3.98 

10.93 

17  Dec. 

2003-17 

March  2004 

28.33 

13.85 

11.34 

8.19 

Izmir  / 
Turkey 

Solar 

assisted 

4.15 

-0.69 

7  Jan.  2004 

23.06 

10.45 

6.53 

9.79 

vertical  GSHP/ 
heating 

(for  greenhouse 
heating) 


21.14 

(127) 

2.94  51.8  23.6  35.4  40.2  [95] 

(129) 

3.84 

17.67 

(127) 

9.07 

[93] 

(17.21) 

(128) 

8.4 

(130) 

66.8 

(131) 

66.6 

30.22 

2.52 

1.83  and 

0.63 

(130) 

71.8 

28.3  9.1  13.5  30.4  [94] 

3.09 

(131) 

67.7 

43.11 

2.61 

1.24  and 

1.30 

(130)  76.2 

23.1  5.6  7.4  21.1  [96] 

1.89 

(131)  75.6 
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Table  8  ( continued ) 


Location 

Type  of 

Capacity 

Dead 

Date 

Exergy  destructions  {rate  of  improvement  potential )  (%  of  the  whole  heat  pump  system) 

Eq.  no.  Exergy 

Thermodynamic 

of  the 

heat  pump/ 

in  the 

state 

of  data 

for  efficiency 

parameters  for 

system 

operating 

mode 

temp. 

used 

Ground  source  heat  pump  system  (Whole  system)3 

exergy  (%) 

compressor  (% 

built/ 

mode 

operated 

(°C) 

eff. 

of  the  whole  system) 

country 

(kW) 

Ground  source  heat  pump  unit  Other  components  Heat  Overall  x  <5  £  / 


-  -  pump 

I  Com.  II  Con.  Ill  Thr.  IV  V  Fan-coil/  VI  GHE  VII  Cir.  VIII 

valve  Evap.  Con.  WT  /PHEx  pumps  Sol. 

col. 


Izmir/ 

Turkey 

Solar  assisted 
vertical  GSHP/ 
heating  (for 
greenhouse 
heating) 

3.97 

1.00  17  Dec. 

2003-17 

March  2004 

3.13 

17.70 

13.68 

6.43 

49.08 

8.61 

0.8 

0.57 

(130) 

(131) 

91.80 

86.13 

[103] 

Erzurum/ 

Turkey 

GSHP  with 
water  source/ 
heating  (for 
space  heating)b 

7.02 

25 

45.92 

15.83 

6.27 

11.69 

14.64 

/5.66 

(130) 

(131) 

4.64 

1.44 

[97] 

aI-Com:  compressor,  II-Con:  condenser,  III-Thr.  valve:  throttle  (expansion)  valve,  IV-Evap:  evaporator,  V:  Fan-coil, /Con.  WT:  condenser  water  tank,  VI-GHE: 
ground-heat  exchanger/PHEx:  plate  heat  exchanger,  VH-Cir.  Pumps:  circulating  pumps,  VIII-Sol.  col.:  solar  collector. 

bThis  system  was  designed  and  constructed  by  Kara  [106]  in  1999  on  the  base  of  a  PhD  thesis  in  the  Mechanical  Engineering  Department,  Ataturk  University,  Erzurum 
in  Turkey  in  order  to  evaluate  geothermal  resources  with  low  temperatures.  The  performance  of  the  system  was  evaluated  by  Kara  and  Yuksel  [107]  using  the  energy 
analysis  method  only,  while  this  system  was  evaluated  by  Akpinar  and  Hepbasli  [97]  using  exergy  analysis  method. 
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A  schematic  diagram  of  the  constructed  experimental  system  is  illustrated  in  Fig.  23 
[11,14].  This  system  mainly  consists  of  three  separate  circuits  as  follows:®  the  ground 
coupling  circuit  with  solar  collector  (brine  circuit  or  water-antifreeze  solution  circuit),  (ii) 
the  refrigerant  circuit  (or  a  reversible  vapor  compression  cycle)  and  (iii)  the  fan-coil  circuit 
for  greenhouse  heating  (water  circuit). 

The  exergetic  results  of  these  two  studies  are  illustrated  in  Table  8.  The  highest 
irreversibility  on  a  system  basis  occurred  in  the  greenhouse  fan-coil  unit,  followed  by  the 
compressor,  condenser,  expansion  valve  and  evaporator,  sub-regions  I  and  V  for  the 
GSHP  unit  and  the  whole  system,  respectively.  Besides  this,  the  remaining  system 
components  have  a  relatively  low  influence  on  the  overall  efficiency  of  the  whole  system. 
Experimental  have  also  showed  that  monovalent  central  heating  operation  (independent  of 
any  other  heating  system)  could  not  meet  overall  heat  loss  of  greenhouse  if  ambient 
temperature  was  very  low.  The  bivalent  operation  (combined  with  other  heating  system) 
could  be  suggested  as  the  best  solution  in  Mediterranean  and  Aegean  regions  of  Turkey,  if 
peak  load  heating  could  be  easily  controlled  [11]. 

4. 3. 2. 4.  Drying.  Large  quantities  of  food  products  are  dried  to  improve  shelf  life, 
reduce  packing  costs,  lower  shipping  weights,  enhance  appearance,  encapsulate  original 
flavor  and  maintain  nutritional  value.  In  evaluating  the  performance  of  food  systems, 
energy  analysis  method  has  been  widely  used,  while  studies  on  exergy  analysis,  especially 
on  exergetic  assessment  of  drying  process,  are  relatively  few  in  numbers  [108-112].  In 
these  previous  studies,  the  drying  process  was  thermodynamically  modeled  by  Dincer 
and  Sahin  [109],  while  drying  of  different  products  such  as  pistachio  [108],  red  pepper 
slices  [110],  potato  [111],  pumpkin  [112],  wheat  kernel  [113]  and  apple  slices  [114]  was 
evaluated  in  terms  of  energetic  and  exergetic  aspects  using  various  drying  devices,  such  as  a 
fluidized  bed  dryer,  a  solar  drying  cabinet,  and  cyclone  type  dryers  [115].  Based  on  the 
previous  studies,  the  first  study  on  exergetic  evaluation  of  a  drying  process  using  a  GSHP 
drying  system  appears  to  be  done  by  Kuzgunkaya  and  Hepbasli  [115],  who  used  laurel 
leaves  as  a  product  being  dried,  and  made  a  comparison  of  exergetic  efficiency  values 
obtained  from  various  studies  on  food  drying  systems.  The  main  objective  of  this 
contribution  was  to  perform  exergy  analyses  of  the  drying  process  in  terms  of  drying  of 
laurel  leaves  under  different  operating  conditions  for  the  assessment  of  the  drying 
performance. 

Fig.  24  [115]  illustrates  a  schematic  diagram  of  this  system,  which  consists  of  mainly 
three  separate  circuits,  namely:  (I)  the  ground  coupling  circuit  (brine  circuit  or 
water-antifreeze  solution  circuit),  (II)  the  refrigerant  circuit  (or  a  reversible  vapor 
compression  cycle)  and  (III)  the  drying  cabinet  circuit  (air  circuit).  The  main  components 
of  the  heat  pump  system  are  an  evaporator,  a  condenser,  a  compressor  and  an  expansion 
valve.  To  avoid  freezing  the  water  under  the  working  condition  and  during  the  winter,  a 
10%  ethyl  glycol  mixture  by  weight  was  prepared.  The  refrigerant  circuit  was  built  on  the 
closed  loop  copper  tubing.  The  working  fluid  was  R-22  [115]. 

There  are  mainly  two  ways  of  formulating  exergetic  efficiency  for  drying  systems 
[108,113,115].  The  first  one  can  be  defined  as  the  ratio  of  the  product  exergy  to  exergy 
inflow  as  follows  [108,115]: 

Exergy  inflow  —  Exergy  loss  Exioss 

^ex  Exergy  inflow  Exin 


(132) 
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Fig.  23.  A  schematic  diagram  of  a  solar  asissted  vertical  ground  source  heat  pump  system  for  greenhouse  heating 
[11,14]. 


Fig.  24.  Schematic  diagram  of  the  heat  pump  tray  drying  system  [116].  (1)  Refrigerant  inlet  to  the  expansion 
valve,  (2)  refrigerant  outlet  from  the  compressor,  (3)  water  inlet  to  the  ground  heat  exchanger,  (4)  water  outlet 
from  ground-heat  exchanger,  (5)  air  inlet  to  the  condenser,  (6)  air  outlet  from  the  condenser,  (7)  air  outlet  from 
the  drying  chamber. 
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The  second  one  may  be  defined  on  the  product/fuel  basis,  as  given  in  Eq.  (133).  The 
product  is  the  rate  of  exergy  evaporation  (Exevap)  and  the  fuel  is  the  rate  of  exergy  drying 
air  entering  the  dryer  chamber  (Exda).  In  this  regard,  exergy  efficiency  may  be  written  as 
follows  [113]. 


Ex 


hex 


evap 


Exd 


a 


with 


Fx  — 
J-^Aevap  — 


l  - 


T 


0 


a 


rp  i  ^evap 

1  m2 


and 


(133) 


(134) 


Sevap  ^w^fg,  (135) 

where  2evap  is  the  heat  transfer  rate  due  to  phase  change,  Tmi  is  the  exit  temperature  of  the 
material,  mw  is  the  mass  flow  rate  of  water  and  Afg  is  the  vaporization  latent  heat. 

The  laurel  leaves  were  sufficiently  dried  at  the  temperatures  ranging  from  40  to  50  °C 
with  relative  humidities  varying  from  16%  to  19%  and  a  drying  air  velocity  of  0.5 m/s 
during  the  drying  period  of  9  h.  The  exergy  efficiency  values  were  obtained  to  range  from 
81.35%  to  87.48%  using  Eq.  (132)  based  on  the  inflow,  outflow  and  loss  of  exergy,  and 
9.11%  to  15.48%  using  Eq.  (133)  based  on  the  product/fuel  basis  between  the  same  drying 
air  temperatures  with  a  drying  air  mass  flow  rate  of  0.12  kg/s  [116]. 

A  comparison  of  exergy  efficiency  values  in  drying  various  products  is  made  in  Table  9 
[108,112,113,115].  It  is  obvious  from  this  table  that  the  exergy  efficiency  values  calculated 
using  Eq.  (132)  are  quite  higher  than  those  using  Eq.  (133).  It  may  be  concluded  that  it  is 
more  meaningful,  objective  and  useful  to  assess  the  performance  of  the  drying  process 
relative  to  the  performance  of  similar  drying  processes  on  the  product  (or  benefit)/fuel 
basis  [116]. 


4.3.3.  Indirect  utilization  (geothermal  power  plants) 

Commercial  delivery  of  geothermally  generated  electric  power  occurred  in  1914  when  a 
250  kW  unit  at  Larderello  provided  electricity  to  the  nearby  cities  of  Volterra  and 
Pomarance  [116].  Since  then,  many  geothermal  power  plants  have  been  installed.  The 
performance  of  these  systems  has  been  mostly  evaluated  using  energy  analysis  method, 
while  the  number  of  the  geothermal  plants,  of  which  performance  has  been  assessed  using 
exergy  analysis  method,  is  very  low.  The  concept  of  exergy  was  first  used  to  analyze  a 
geothermal  power  plant  by  Badvarsson  and  Eggers  [117].  Their  illustrative  example,  which 
compared  the  performances  of  single  and  double  flash  cycles  based  on  a  reservoir  water 
temperature  of  250  °C  and  a  sink  condition  of  40  °C,  gave  exergetic  efficiencies  of  38.7% 
and  49%,  respectively,  assuming  65%  mechanical  efficiency  [71]. 

Using  Eq.  (15)  along  with  its  concept  and  applying  it  to  a  power  plant  as  a  whole,  the 
overall  exergetic  efficiency  reduces  to  a  very  simple  formula,  namely,  the  ratio  of  the  net 
power  output  ( JUnet, plant)  to  the  exergy  of  the  motive  fluid  serving  as  the  energy  source  for 
the  plant  [118].  That  is: 

_  E  net, plant 

%pp  “  Exwh 


(136) 


Table  9 

Comparison  of  exergy  efficiency  values  calculated  by  some  investigators  in  drying  various  products  [108,112,113,115] 


Investigators/year  of 
publication 

Type  of  dryer 

Product  dried 

Air 

Exergy  efficiency 

Type 

Initial 
moisture 
content  (%) 

Final 

moisture 

content  on 
dry  weight 
basis  (%) 

Velocity  (m/s) 

Temperature 

(°C) 

Relative 
humidity  (%) 

Eq.  no. 

Value  (%) 

Syahrul  et  al.  [113] 

Fluidized  bed 

Red-spring 

31.7-32.6 

15 

1.95 

40.2-65 

18.5-21.1 

(133) 

4-12 

dryer 

wheat 

Corn 

24.6-25.6 

2.22-2.24 

50-63 

15.2-17.5 

2-16 

Midilli  and  Kucuk 

Solar  drying 

Shelled/ 

26.95/29 

1.23 

40-60 

37-62 

(132) 

10.86-100 

[108] 

cabinet 

unshelled 

pistachios 

Akpinar  et.  al.  [112] 

Cyclone  type 

Pumpkin 

6 

1.5 

60 

10-20 

(132) 

30.81-100 

dryer 

slices 

6 

70 

37.89-100 

4 

80 

46.97-100 

Kuzgunkaya  and 

Ground- 

Laurel  leaves 

94.4 

12 

0.5 

40 

16-19 

(132) 

81.35-87.48 

Hepbasli  [115] 

source  heat 

50 

(133) 

9.11-15.48 

pump  dryer 

40 

50 
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where  Exwh  and  Exres  are  the  total  exergy  inputs  to  the  plant  at  wellhead  and  reservoir 
conditions,  respectively,  while  they  are  calculated  using  Eqs.  (6c)  and  (7). 

Note  that  the  exergy  input  to  the  plant  is  sometimes  taken  as  the  exergy  of  geothermal 
fluid  in  the  reservoir.  Those  who  prefer  this  approach  argue  that  a  realistic  and  most 
meaningful  comparison  between  geothermal  power  plants  can  only  be  performed 
considering  the  methods  of  harvesting  the  geothermal  fluid.  However,  others  argue  that 
taking  the  reservoir  as  the  starting  point  is  not  proper  for  geothermal  power  plants  since 
conventional  power  plants  are  evaluated  on  the  basis  of  the  exergy  of  the  fuel  burned  at  the 
plant  site  [118,119]. 

Exergetic  efficiencies  of  some  geothermal  power  plants  worldwide  in  order  of  increasing 
efficiency  are  summarized  by  the  author,  as  illustrated  in  Table  10  [119,120-127].  It  is 
obvious  from  this  table  that  exergy  efficiencies  vary  from  16.3%  to  53.9%  depending  on 
the  dead  state  temperature  and  the  technology  used. 

DiPippo  [120]  obtained  the  following  correlation  for  the  exergetic  efficiency  as  a 
function  of  the  dead-state  temperature  using  the  Brady  geothermal  power  plant: 


eBrady  =  23.177  -  0.14297%  (138) 

DiPippo  [120]  concluded  as  follows:  “Geothermal  binary  plants  are  relatively  poor 
converters  of  heat  into  work.  Thermal  efficiencies  typically  lie  in  the  range  of  8-12%.  As  a 
consequence,  a  1-2%  point  improvement  in  power  output  translates  into  a  gain  of  about 
10-20%  in  efficiency.  The  results  show  that  binary  plants  can  operate  with  very  high 
exergetic  efficiencies  even  when  the  motive  fluids  are  low-temperature  and  low-exergy. 
Exergetic  efficiencies  of  40%  or  greater  have  been  achieved  in  certain  plants  with  geofluids 
having  specific  exergies  of  200kJ/kg  or  lower.  The  main  design  feature  leading  to  a  high- 
exergy  efficiency  lies  in  the  design  of  the  heat  exchangers  to  minimize  the  loss  of  exergy 
during  heat  transfer  processes.  Another  important  feature  that  can  result  in  a  high  exergy 
efficiency  is  the  availability  of  low-temperature  cooling  water  that  allows  a  once-through 
system  for  waste  heat  rejection.” 

In  the  following,  the  Grassmann  diagram  of  a  12.4  MW  existing  binary  geothermal 
power  plant  is  presented  as  an  illustrative  example  of  geothermal  power  plants.  Exergy 
analysis  of  this  plant  was  performed  by  Kanoglu  [119]  using  actual  plant  data  to  assess  the 
plant  performance  and  the  Grassmann  shown  in  Fig.  25  [119]  was  drawn.  As  seen  in  this 
figure,  the  causes  of  exergy  destruction  in  the  plant  included  the  exergy  of  the  working  fluid 
lost  in  the  condenser,  the  exergy  of  the  brine  reinjected,  the  turbine -pump  losses,  and  the 
preheater-vaporizer  losses.  The  exergy  destruction  at  these  sites  accounted  for  22.6%, 
14.8%,  13.9%,  and  13.0%  of  the  total  exergy  input  to  the  plant,  respectively. 

Dagdas  et  al.  [127]  proposed  a  new  model  in  order  to  improve  efficiency  and  augment 
power  output  of  the  single  flash  type  geothermal  power  plant  operating  in  Denizli- 
Kizildere,  Turkey.  In  this  model,  electric  energy  requirement  of  the  region  would  be  met 
and  performance  of  the  geothermal  power  plant  would  be  up  rated  by  building  a  gas 
turbine  power  plant  near  the  existing  plant.  If  the  proposed  combined  cycle  power  plant 
(gas  turbine  integrated  geothermal  power  plant)  would  be  operated,  power  generated  in  the 
geothermal  plant  could  increase  43.5%.  Energy  and  exergy  efficiencies  of  the  plant  would 
be  increased  from  4.56%  to  6.5%,  and  from  19.78%  to  28.4%,  respectively. 


Table  10 

Geothermal  power  plant  exergetic  efficiencies  (in  order  of  increasing  efficiency)  [119,120-127] 


Technology 

Plant  name 

Location/country 

Dead  state  temperature 

(°C) 

Exergy  efficiency a  (%) 

Sources 

Binary 

Brady  bottoming 
(p.m.  data) 

Churchill  County, 
NV/USA 

16.8 

16.3 

DiPippo  [120] 

Binary 

Brady  bottoming 
(a.m.  data) 

Churchill  County, 
NV/USA 

30.1 

17.9 

DiPippo  [120] 

Binary:  recuperated 

Rotokawa 

Rotokawa/New  Zealand 

18.7 

DiPippo  [120] 

Single  flash 

Kizildere 

Denizli/Turkey 

16 

19.78a 

19.97b 

Dagdas  et  al.  [121,127] 

Single  flash 

Kizildere 

Denizli/Turkey 

16 

20.5 

Ozturk  et  al.  [122] 

Single  flash 

Kizildere 

Denizli/Turkey 

15 

20.8 

Cerci  [123] 

Single  flash 

Kizildere 

Denizli/Turkey 

18 

24. 3b 

27. 2C 

Yildirim  and  Gokcen 
[124] 

Binary 

Nigorikawa  (Mori)  pilot 

Hakodate  on 
Hokkaido/Japan 

13 

21.6 

DiPippo  [120] 

Binary 

Kalina  KCS  34,  Husavik 

Husavik/  Iceland 

5 

23.1 

DiPippo  [120] 

Single  flash 

Kakkonda 

Japan 

26.5 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

Binary:  simple 

Rotokawa 

Rotokawa/New  Zealand 

27.8 

DiPippo  [120] 

Single  flash 

Otake 

Japan 

29.0 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

Dual-level  binary 

Northern  NV/USA 

12.8 

29. ld 

34.2e 

Kanoglu  [119] 

Single  flash 

Salak  (Unit  1&2) 

Indonesia 

30.4 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 
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F 

Binary:  simple 

Kavala 

Eratino/Greece 

Single  flash 

Ahuachapan 
(Unit  1  &  2) 

El  Salvador 

Single  and  double  flash 

Cerro-Prieto 

Mexico 

Single  and  double  flash 

Ahuachapan 
(Unit  1  &  2) 

El  Salvador 

Single  flash 

Blundell 

Near  Milford  in  Beaver 
County/USA 

Hybrid  flash-binary 

Rotokawa 

Rotokawa/New  Zealand 

Binary:  dual-level 

Heber  SIGC 

Imperial  Valley,  CA/ 
USA 

Double  flash 

Beowawe 

NV/USA 

Dry  steam 

Kamojang 
(Unit  2  &  3) 

Indonesia 

Dry  steam 

Darajat 

Indonesia 

Binary:  flash  evaporator 

Otake  pilot 

Island  of  Kyushu/ Japan 

20 

33. 3g 

8.1g 

52. 2h 

Koroneos  et  al.  [126] 

33.1 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

33.2 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

38.2 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

35.6 

DiPippo  [120] 

42.0 

DiPippo  [120] 

15 

43.4 

DiPippo  [120] 

46.0 

DiPippo  [120] 

48.4 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

52.4 

Setiawan  [125],  Yildirim 
and  Gokcen  [124] 

18 

53.9 

DiPippo  [120] 

aExergy  efficiencies  except  those  denoted  in  the  above  table  are  based  on  reservoir  exergy. 
bAt  reservoir  conditions. 
cAt  wellhead  conditions. 

d Based  on  the  exergy  of  the  geothermal  fluid  at  the  vaporizer  inlet. 

eBased  on  the  exergy  drop  of  the  brine  across  the  vaporizer-preheater  system  (i.e.,  exergy  input  to  the  Rankine  cycle). 

f 

For  the  exploitation  of  the  field,  a  Rankine  cycle  is  proposed  by  Koroneos  et  al.  [126]. 
gFor  the  first  deposit  based  on  the  two  case  studies. 
hFor  the  second  deposit. 
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Brine 

reinjection 


Fig.  25.  Grassmann  (exergy  flow)  diagram  of  a  12.4  MW  existing  binary  geothermal  power  plant  drawn  by 
Kanoglu  [119]. 


4.4.  Biomass 


Chemical  exergy  of  a  substance  is  the  maximal  possible  useful  work  that  may  be 
produced  by  process  of  physical  and  chemical  equilibration  of  the  substance  with  the 
ambient.  To  evaluate  the  chemical  exergy  of  a  mixture,  the  knowledge  of  its  enthalpy  of 
combustion,  elemental  composition  and  absolute  entropy  is  necessary.  When  some  of  these 
data  are  absent,  the  methods  of  estimation  of  energy  are  used  [128].  The  chemical  exergy  of 
gaseous  and  some  liquids  can  be  obtained  using  standard  chemical  exergy  tables  given  in 
Ref.  [7]  if  the  chemical  composition  of  the  fuel  is  known.  Many  technical  liquid  and  solid 
fuels  are  multi-component  solutions  (or  mixtures)  of  very  complicated,  usually  unknown 
components.  An  exact  solution  of  the  chemical  exergy  of  such  fuels  is  impossible.  In  this 
regard,  Rant’s  relation  including  a  constant  ratio  of  chemical  exergy  and  calorific  value  (or 
heating  value)  for  liquid  and  (separately)  for  solid  fuels  (Exch/HV)  may  be  used. 
However,  the  calculation  for  different  organic  substances  indicated  that  the  ratio 
(Exch/HV)  depends  significantly  on  the  chemical  composition. 

If  the  elemental  composition  is  known,  the  chemical  exergy  of  various  fuels  can  be 
evaluated  accurately  using  the  Szargut’s  correlations  listed  elsewhere  [7].  In  addition,  a 
correlation  formula  to  estimate  chemical  exergies  of  oil  fractions  and  fuel  mixtures  from 
enthalpy  of  combustion  and  atomic  composition  developed  by  Govin  et  al.  [128]  is  used. 

Based  approximation  method  for  fuel  chemical  exergy  (<Ach),  the  following  relations  are 
widely  used. 


/^LHV  — 


ExCh 


LHV  ’ 


(139) 
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^hhv~HHV’  (140) 

where  /i  is  the  proportionality  constant  (or  quality  factor  or  exergy  coefficient),  while  LHV 
and  HHV  denote  the  lower  heating  (net  caloric)  and  higher  heating  (gross  caloric)  values, 
respectively  [129]. 

Szargut  [22]  proposed  that  /?LHv  and  /?hhv  may  be  taken  to  be  approximately  1.15  and 
1.05  for  wood,  respectively.  The  following  relations  for  wood  with  the  limit  of 
(Z02/ZcK2  .67  are  proposed  by  Szargut  et  al.  [7]  and  Szargut  [22],  respectively. 

£lhv  =  {1.0412  +  0.2160(Zh,/Zc)  -  0.2499(ZO2/Zc)[l  +  0.7884(ZH2/ZC)] 

+0.0450(Zn2/Zc)}/{  1  -  0.3035(Zo2/Zc)}  (141) 

and 

£lhv  =  {1.0412  +  0.2160(Zh2/Zc)  +  0.2499(ZO2/Zc)[l  +  0.7884(ZH2/ZC)] 

+0.0450(Zn2/Zc)}/{  1  +  0.3035(Zo2/Zc)},  (142) 

where  Z?  is  mass  fraction  of  element  i. 

Eq.  (141)  has  been  successfully  used  by  some  investigators  [130-132].  In  this  regard, 
Nilsson  [130]  used  the  following  relation  for  determining  the  specific  chemical  exergy  of 
straw 


E^CH, straw  —  (LHV  +  TwZw)/?m y  T  Exch,w^w?  (143) 

where  L  and  Z  denote  the  enthalpy  of  phase  change  (kJ/kg)  and  the  mass  fractions  ( — ), 
respectively.  Using  the  following  data  for  moist  straw;  Zq  =  38.9%,  Zh,  =  4.9%, 
Z02  =  33.8%,  ZN,  =  0.4%,  Zw  =  18.0%,  LHV  =  13,800  kJ/kg,  Lw  =  2440 kJ/kg  and 
ExCH,w=  50  kJ/kg,  the  specific  chemical  exergy  was  calculated  to  be  16,100kJkg.  In 
other  words,  the  ratio  of  chemical  exergy  to  the  lower  heating  value  was  1.16.  In  the 
calculation,  the  reference  temperature  was  25  °C,  the  relative  humidity  of  atmospheric  air 
was  70%  and  the  chemical  exergy  of  ash  was  neglected. 

Panopoulos  et  al.  [132]  analyzed  exergetically  the  combination  of  allothermal  biomass 
gasification  and  solid  oxide  fuel  cell  for  small  scale  CHP,  and  presented  exergy  efficiency 
relations.  They  also  determined  that  for  SOLC  operation  at  current  density  2500  A/m  ,  the 
system  used  90kg/h  biomass  and  operated  with  electrical  exergetic  efficiency  32% 
producing  140  kWe,  while  the  combined  electrical  and  thermal  exergetic  efficiency  was 
35%. 

Lranco  and  Giannini  [133]  analyzed  the  perspective  of  improving  the  efficiency  of 
biomass  for  power  generation  and  proposed  an  exergy  loss  based  economic  analysis  of 
biomass  utilization.  A  thermoeconomic  method  of  analysis  based  on  the  consideration  of 
the  cost  of  the  exergy  losses  showed  how  the  conversion  efficiency  of  renewable  sources  as 
biomass,  had  a  lower  limit  that,  in  a  lot  of  cases,  was  considerably  higher  than  the  actual 
efficiency  level. 

According  to  Shieh  and  Lan  [134],  specific  chemical  exergy  (kJ/kg)  of  structurally 
complicated  materials,  e.g.  biomass,  can  be  estimated  from  their  elemental  compositions  as 
given  below  [135] 

t/,CH  =  4.1868{8177.79[C]  +  5.25[N]  +  27,892.63[H]  -  3173.66(0] 

+  0. 1 5[0](7837.677[C]  +  33, 888.889[H]  -  4236.1[0]}. 


(144) 
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Table  11 

The  calculated  exergy  values  of  various  of  woody  biomass  types  in  Turkey  [136] 


Type  of  wood  Ultimate  analysis  (%)  Exergy  (MJ/kg)  HHV  (MJ/kg)  /?HHV 


CHON 


Hus  wood 

48.5 

5.6 

45.7 

0.2 

17.90 

17.22 

1.040 

Willow  wood 

50.8 

6.0 

42.5 

0.6 

19.63 

19.01 

1.033 

Sugar  cane 

47.3 

6.2 

46.2 

0.3 

18.15 

17.52 

1.036 

Olive  residue 

49.5 

6.3 

43.7 

0.5 

19.39 

18.77 

1.033 

Straw 

47.0 

6.1 

46.4 

0.5 

17.90 

17.25 

1.038 

Corn  cob 

46.2 

6.3 

47.2 

0.3 

17.75 

17.12 

1.037 

Olive  shell 

52.4 

6.0 

40.4 

1.5 

20.45 

19.76 

1.031 

Wheat  straw 

48.7 

6.2 

44.8 

0.3 

18.83 

18.21 

1.034 

Rice  shall 

48.7 

6.0 

44.7 

0.6 

18.60 

17.96 

1.036 

Rice  straw 

47.0 

6.4 

45.5 

1.1 

18.38 

17.77 

1.034 

Hybride  poplar 

51.6 

6.3 

41.5 

0.6 

20.42 

19.84 

1.029 

Cotton  stalk 

46.7 

6.2 

46.5 

0.6 

17.90 

17.96 

1.037 

Bilgen  et  al.  [136]  studied  on  exergy  analyses  of  various  types  of  wood  and  coal  samples 
with  Turkish  origins  by  using  exergy  analysis  method  developed  and  published  by 
Srivastava  [137]  and  Bejan  et  al.  [138].  They  calculated  the  exergy  values  of  various  of 
woody  biomass  types  in  Turkey,  as  given  in  Table  11  [136].  As  can  be  seen  in  this  table,  the 
values  for  /?Hhv  vary  from  1.029  to  1.040.  On  the  other  hand,  these  values  were  found 
greater  than  1.0  for  all  the  woody  biomass  types.  Bilgen  et  al.  [136]  concluded  that  this 
result  proposed  that  maximum  work  is  always  greater  than  the  heating  value,  or  energy 
delivered  by  a  process  in  which  these  coals  are  used  as  a  functional  material. 

Utlu  et  al.  [139]  evaluated  the  Turkish  cotton  stalk  production  in  2003  using  energy  and 
exergy  analysis  method  for  the  first  time.  The  values  for  the  net  energy  and  exergy  gained 
were  obtained  to  be  about  49,146  and  38,426  MJ/ha,  respectively.  The  overall  mean  energy 
and  exergy  efficiencies  of  the  cotton  production  in  the  year  studied  were  found  to  be 
33.06%  and  33.12%,  respectively. 

Prins  et  al.  [140]  performed  laboratory  experiments  in  a  small-scale  (5-10  g  sample) 
fixed-bed  torrefaction  reactor  placed  inside  an  oven  and  gave  exergetic  efficiency 
relations  for  various  processes.  They  reported  as  follows:  Overall  exergetic  efficiency  of 
air-blown  gasification  of  torrefied  wood  was  lower  than  that  of  wood  since  the  volatiles 
produced  in  the  torrefaction  step  were  not  utilized.  For  the  entrained  flow  gasifier,  the 
volatiles  could  be  introduced  into  the  hot  product  gas  stream  as  a  ‘chemical  quench’.  The 
overall  efficiency  of  such  a  process  scheme  was  comparable  to  direct  gasification  of  wood, 
but  more  exergy  was  conserved  in  as  chemical  exergy  in  the  product  gas  (72.6%  versus 
68.6%). 


4.5.  Other  renewable  energy  systems 

The  following  section  summarizes  exergy  relations  for  ocean  surface  waves,  pre¬ 
cipitation,  ocean  thermal  gradient  and  tides,  which  have  been  taken  from  a  study 
conducted  by  Hermann  [14].  These  relations  may  be  used  in  exergetic  evaluating  these 
types  of  RERs. 
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4.5.1.  Ocean  surface  waves 

The  average  specific  exergy  of  surface  waves  ranges  from  10  to  lOOkW/m,  with  much 
higher  values  during  storms  [14,141].  The  waves  concentrate  kinetic  energy  transfer  from 
the  wind  over  many  kilometers  and  can  deposit  this  energy  at  the  coastline.  The  primary 
factors  determining  the  exergy  of  a  fluid  gravity  wave  are  the  period  and  [14,142].  For  an 
approximately  sinusoidal  deep  ocean  wave,  the  length-specific  exergy  flux  is 

1  7 

*Awa  ve  2^^wave^wave’  (145) 

where  Fwave  is  the  wave  speed,  g  is  the  gravitational  constant  and  aWSLWG  is  the  amplitude. 


4.5.2.  Precipitation 

The  specific  exergy  of  precipitation  can  be  broken  down  into  physical  and  chemical 
components  due  to  its  gravitational  and  diffusive  chemical  potential  energy  relative  to 
seawater. 

For  a  specific  replenishing  body  of  water  such  as  a  dammed  reservoir,  the  gravitational 
potential  exergy  is 

'/'pre.pot  =  (146) 

where  m  is  the  mass  flow  rate  and  z  is  the  vertical  distance  from  the  water  level  of  the 
reservoir  to  the  reference  height  or  average  sea  level  [14]. 

Precipitation  also  contains  a  diffusive  exergy  component  due  to  the  difference  in 
chemical  potential  between  the  Earth’s  freshwater  and  seawater.  Since  the  difference 
between  seawater  and  freshwater  is  primarily  dissolved  salts,  this  value  can  be  computed 
most  directly  by  considering  a  reversible  osmosis  process. 

The  work  done  against  the  osmotic  pressure  is 

•Apre.diff  =  -RT  ln(l  -  .y solute)’  (147) 

where  y  is  the  mole  fraction  of  solute  in  the  water  [14,143]. 


4.5.3.  Ocean  thermal  gradient 

The  solar  warming  of  seawater  in  the  tropics  and  the  radiative  loss  of  thermal  energy  to 
space  near  the  poles  creates  a  density-driven  current  which  circulates  seawater  from  the 
surface  to  the  deep  ocean.  Due  to  the  small  temperature  differential,  the  exergy  density  of 
the  ocean  thermal  gradient  is  relatively  small.  For  a  temperature  difference  of  20  K,  the 
exergy  content  is  approximately  7%  of  the  thermal  energy  [14]. 


*AtEG  —  ^fi^p,ave(TWarm  —  Tcoid)  (  1  —  —  j  . 

V  1  warm  / 


(148) 


4.5.4.  Tides 

The  specific  tidal  exergy  of  a  designated  area  of  ocean  or  enclosed  reservoir  is  a  function 
of  the  record  of  minima  and  maxima  and  the  time  periods  between  these  levels  [14].  As  the 
tide  rises  and  falls,  an  average  exergy  flow  considering  the  reservoir  area  integrated  over 
the  height  difference  between  local  maxima  and  minima  and  the  period  between  these 
points  can  be  written  as  follows: 
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where  t  is  the  period  between  local  maxima  and  minima  of  the  tidal  record,  Ares  the 
area  of  the  reservoir  at  height  z,  and  n  the  index  number  for  each  maximum  or 
minimum  [14]. 


4.6.  Country  based  renewable  energy  sources 

The  method  of  exergy  analysis  has  been  applied  to  a  wide  variety  of  thermal  and  thermo 
chemical  systems.  A  particular  thermo  dynamical  system  is  the  society,  for  example,  of  a 
country  or  a  region  [144].  In  order  to  attain  energy  saving,  the  application  of  energy  and 
exergy  analysis  modeling  techniques  for  energy-utilization  assessments  has  recently 
significantly  increased.  Utlu  and  Hepbasli  [145]  conducted  a  comprehensive  review  on  the 
exergetic  assessment  of  countries.  In  this  regard,  the  approaches  used  to  analyze  energy 
utilization  of  countries  or  societies  may  be  grouped  into  three  types,  namely  Reistad’s 
approach.  Waifs  approach  and  Sciubba’s  approach.  In  Reistad’s  approach,  the  end  use  is 
divided  into  three  sectors,  these  names  are  industry,  transportation  and  residential-com¬ 
mercial.  In  addition,  energy  sector  with  electric  generation  and  distribution,  oil  refining,  is 
treated  separately.  Flows  of  energy  carriers  for  energy  uses  are  not  included  in  this 
approach.  In  Wall’s  approach,  all  types  of  energy  and  material  flows  are  considered. 
However,  to  energy  carriers  for  energy  use  are  considered,  these  flows  encompass  wood  for 
construction  materials  and  for  the  pulp  and  paper  industry,  harvested  wood  and  fodder, 
and  the  products  from  these  materials.  The  so-called  Sciubba’s  approach  is  relatively  a  new 
one.  Extended-exergy  accounting  (EEA)  method  introduced  by  Sciubba  was  applied  to  the 
Italian  society  in  1996  by  Milia  and  Sciubba  [146]  and  recently  to  Norway  with  the  figures 
of  2000  by  Ertesvag  [147]. 

Table  12  [145,149,150]  illustrates  Eqs.  (1 50)— (1 58)  used  to  estimate  exergy  (second  law) 
utilization  efficiencies  for  the  RERs.  The  exergy  efficiency  frequently  gives  a  finer 
understanding  of  performance  than  (g).  The  parameter  (g)  weights  energy  flows  by 
accounting  for  each  in  terms  of  availability.  It  stresses  that  both  losses  and  internal 
irreversibilities  need  to  be  dealt  with  to  improve  performance. 

To  date,  various  studies,  as  reviewed  in  detail  elsewhere  [145],  were  performed  to 
analyze  energy  and  exergy  utilization  efficiencies  of  countries  in  the  four  subsectors, 
namely  utility,  industrial,  residential-  commercial  and  transportation  sectors.  However, 
studies  conducted  on  the  evaluation  of  energy  utilization  efficiency  of  RERs  using  energy 
and  exergy  analyses  are  very  limited  in  numbers  in  the  open  literature.  The  works  done  by 
Hepbasli  and  Utlu  [15,148]  on  Turkey’s  RERs  for  the  years  of  2001  and  2010  appear  to  be 
the  first  ones. 

Table  13  [15,148,151-154]  illustrates  the  main  outcomes  of  the  exergy  utilization  analysis 
for  both  Turkey’s  general  and  RERs.  It  is  clear  from  this  table  that  the  exergy  utilization 
efficiencies  for  Turkey’s  general  ranged  from  13.10  to  27.94%.  For  Turkey’s  RERs,  exergy 
utilization  efficiencies  varied  from  24.14  and  44.30%  in  2001  and  2010,  respectively  [15]. 
These  results  have  clearly  indicated  the  necessity  of  the  planned  studies  towards  increasing 
exergy  utilization  efficiencies  in  the  sector  studied.  Hepbasli  and  Utlu  [15]  concluded  that  the 
analyses  reported  here  would  provide  the  investigators  with  knowledge  about  how  much 
effective  and  efficient  a  country  uses  its  natural  resources,  while  this  knowledge  is  also  needed 
for  identifying  energy  efficiency  and/or  energy  conservation  opportunities  as  well  as  for 
dictating  the  energy  strategies  of  a  country  or  a  society,  especially  for  a  sustainable  future. 


Table  12 

Relations  used  in  the  exergy  analysis  [145,149,150] 


Equations  Description 


1  £  =  (exergy  in  products/total  exergy  input)  x  100  (150) 


2  £  =  (r,/P){  1  -  [T0/(T-  T0)\\n(T/T0)}  (151) 

3  £  =  (^out/^in)100  (152) 

4  s  =  rj[\-(T0/T)]  (153) 

5  s  =  ri[(T0/T)-\]  (154) 

6  £oe  =  [(«l£al)  +  («2Sa2)  +  (^3^)  H - h  («7£a7)]/100  (155) 

7  fiof  [(fsh  x  £sh,f)  +  (/wh  x  £wh,f)  +  (/ c  X  £c,f)]/100  (156) 

8  ^or  —  [(^: sh  X  £2sh,or)  T  (r  wh  X  £wh,orf)  T  (^c  X  £2c,ro)]/(s  sh  T  £  wh  T  rc)  (157) 

9  So  —  (Se^rc  T  Sot/* exrc  T  Sr^exrc)/(^rc  T ,/'exrc  T  ^exrc)  (158) 


Definition  of  exergy  utilization  efficiencies 

Exergy  utilization  efficiencies  of  space  heating 

Exergy  utilization  efficiencies  of  RERs  for  utility  sector 

Exergy  utilization  efficiencies  of  cooking 

Exergy  utilization  efficiencies  of  refrigeration 

Overall  sector  exergy  utilization  efficiencies  for  electric  use 

Overall  sector  exergy  utilization  efficiencies  for  fuel  use 

Overall  sector  exergy  utilization  efficiencies  for  RE  use 

Overall  sector  exergy  utilization  efficiencies  for  residential-commercial  sector 
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Table  13 

Comparison  of  total  exergy  inputs/outputs  as  well  as  exergy  utilization  efficiencies  of  Turkey’s  renewables  and  total  [15,148,151-154] 


Investigators 

Year 

analyzed/ 

published 

Population 

Total  exergy  input 

Total  exergy  output 

RERsc 

inputa 

RERsc 

output 

Efficiencies 

(pj) 

(GJ/cap) 

(pj) 

(GJ/cap) 

(pj) 

(GJ/cap) 

(pj) 

(GJ/cap) 

Total 
£  (%) 

RERsc 

£  (%) 

Unal  [151] 

1991/1994 

57,024,515 

2279 

39.97 

539.5 

9.46 

450.54 

7.90 

23.7 

Rosen  and  Dincer  [152] 

1993/1997 

58,808,625 

680.6 

11.57 

445.2 

7.57 

455.1 

7.74 

27.1 

Ileri  and  Gurer  [153] 

1995/1998 

59,706,545 

2697.3 

45.17 

352.3 

5.9 

474.4 

7.95 

13.  lb 

Utlu  and  Hepbasli  [154] 

1999/ 

66,022,636 

3380.34 

51.2 

499.61 

7.55 

459.25 

6.96 

24.04 

2000/ 

67,803,927 

3469.62 

52.55 

525.77 

7.96 

449.72 

6.63 

24.78 

Hepbasli  and  Utlu  [15] 

2001/2004 

68,820,985 

3139.07 

45.61 

783.75 

11.39 

395.68 

5.75 

95.51 

1.38 

25.22 

24.14 

Hepbasli  and  Utlu  [148] 

2010/ 

76,840,418 

7400.34 

96.31 

2067.65 

26.91 

580.21 

7.55 

351.7 

4.58 

27.94 

44.30 

aThese  values  are  found  by  calculation  from  references. 

bExcluding  utility  sector  use.  These  values  are  obtained  to  be  43.61  (21.83)%  with  the  utility  sector  use  of  236.6 PJ. 
cRenewable  energy  resources. 
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5.  Conclusions 

Exergetic  aspects  of  RERs  are  comprehensively  presented  in  this  study.  These  RERs 
studied  are  solar,  wind  and  geothermal  energy  systems  as  well  as  biomass  and  country 
based  RERs.  Studies  conducted  on  these  RERs  are  also  compared  with  the  previously  ones 
in  tabulated  forms. 

Some  concluding  remarks  which  can  be  extracted  from  this  study  are  as  follows: 

(a)  Exergy  is  a  way  to  a  sustainable  development.  In  this  regard,  exergy  analysis  is  a  very 

useful  tool,  which  can  be  successfully  used  in  the  performance  evaluation  of  RERs  as 

well  as  all  energy-related  systems. 

(b)  Based  on  the  review  results  of  the  RERs,  the  following  exergy  efficiency  values  were 

found. 

•  The  exergy  efficiency  of  a  SPC  was  found  to  range  from  0.4%  to  1.25%  [33]. 

•  The  exergy  efficiencies  of  a  solar  collector,  a  PV  and  a  hybrid  solar  collector 
were  found  to  be  4.4%,  11.2%  and  13.3%,  respectively.  Solar  collector  had 
the  lowest  value.  In  terms  of  energy  quality,  the  hybrid  solar  collector  was  very 
effective  [67]. 

•  Based  on  an  exergetic  evaluation  of  wind  energy  systems,  average  differences 
between  energy  and  exergy  efficiencies  were  obtained  to  be  approximately  40%  at 
low  wind  speeds  and  up  to  approximately  55%  at  high  wind  speeds  [68]. 

•  Geothermal  resources  should  be  classified  by  their  exergy,  a  measure  of  their  ability 
to  do  work.  The  classification  with  reference  to  SExI  was  more  meaningful  as  there 
was  no  general  agreement  on  the  arbitrary  temperature  ranges  used  in  the 
classification  of  geothermal  resources  by  temperature  [70,71]. 

•  Based  on  a  comparison  of  exergy  efficiency  values  for  GDHSs  installed  in  Turkey, 
exergy  efficiencies  ranged  from  42.94%  to  66.13%  depending  on  the  various  dead 
state  temperatures.  As  expected,  the  lower  the  reference  state  temperature,  the 
significantly  larger  the  exergy  losses  in  the  system  pipeline  and  heat  exchangers 
[81-83,86,87,90,91]. 

•  Exergy  efficiency  values  of  GSHP  units  varied  between  about  67%  and  92%  on 
the  product/fuel  basis  depending  on  the  various  dead  state  temperatures 
[93-97,103,106,107]. 

•  It  was  found  from  an  exergetic  evaluation  of  a  drying  process  utilizing  a  GSHP 
system  that  the  exergy  efficiency  values  of  the  drying  process  ranged  from  81.35%  to 
87.48%  based  on  the  inflow,  outflow  and  loss  of  exergy,  and  9.11  %  to  15.48%  based 
on  the  product/fuel  basis  between  the  same  drying  air  temperatures  with  a  drying  air 
mass  flow  rate  of  0.1 2  kg/s.  It  was  also  concluded  that  it  was  more  meaningful, 
objective  and  useful  to  assess  the  performance  of  the  drying  process  relative  to  the 
performance  of  similar  drying  processes  on  the  product  (or  benefit)/fuel  basis  [115]. 

•  Exergetic  efficiencies  of  some  geothermal  power  plants  worldwide  were  compared 
with  each  other  and  found  to  vary  from  16.3%  to  53.9%  depending  on  the  dead 
state  temperature  and  the  technology  used  [119,120,122-127]. 

•  Exergy  utilization  efficiencies  of  Turkey’s  RERs  were  found  to  range  from  24.14% 
and  44.30%  in  2001  and  2010,  respectively  [15].  These  results  have  clearly  indicated 
the  necessity  of  the  planned  studies  towards  increasing  exergy  utilization  efficiencies 
in  the  sector  studied. 
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(c)  As  a  conclusion,  the  author  expects  that  the  analyses  reported  here  will  provide  the 
investigators  with  knowledge  about  how  effective  and  efficient  a  country  uses  its  RESs. 
This  very  useful  knowledge  is  also  needed  for  identifying  energy  efficiency  and/or 
energy  conservation  opportunities,  as  well  as  for  dictating  the  right  energy  and  exergy 
management  strategies  of  a  country. 
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